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^  A  generalized  analytical  model  a  f  the  wind,  current,  ana  wave  forces  on  an  ellipti¬ 
cal  cylinder  is  presented  with  specific  application  to  predicting  the  list  angle  of  \ 
an  articulated  spar  buoy.  The  major  physical  and  environmental  parameters  present 
in  the  system  are  then  analyzed  to  determine  their  relative  influence  on  these 
forces  and  moments. 

Data  to  validate  this  analytical  mocel  was  obtained  from  both  laboratory  and  field 
experiments,  'todei  tests  performed  in  a  circulating  water  cnannel  revealed  that 
directional  stability  was  a  significant  factor  in  determining  the  overall  perfor¬ 
mance  of  elliptically  cross-sectioned  spars.  Of  the  shapes  investigated,  only  a 
circular  cyl  inder  with  splitter  plate  was  found  to  be  staple  ana  to  return  to  an 
equilibrium  position  if  displaced.  Laboratory  measurements  of  the  circular  cylinder 
list  angle  far  various  current  velocities  showed  generally  good  agreement  with  that 
predicted  by  the  analytical  model. 

Using  a  nominal  6  inch  (16.83  cn)  diameter,  in;.  '  aluminum  spar  18  feet  (5.4S 

m)  in  length  with  a  splitter  plate,  field  "esc  <  i  ootained  over  a  broad  -ange 

of  current  and  water  depth  conditions.  3ata  from  thi.  eld  experiment  also  showed 
good  agreement  with  that  predicted  Sy  the  analytical  model. 

On  the  basis  of  these  results,  it  was  concluded  that  for  soars  of  circular  cross- 
section,  the  analytical  model  is  a  good  predictor  of  an  articulated  soar's  benavior 
and,  consequently,  a  useful  design  aid.  Possible  improvements  to  both  the  laoora- 
tory  and  field  tests  are  suggested. 
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CHAPTER  I 


BACKGROUND 

1-1  Introduction 

The  U.S.  Coast  Guard  presently  has  over  24,000  lighted  and  unlighted 
buoys  serving  as  aids  to  marine  navigation  in  harbors,  rivers,  lakes  and 
coastal  areas  across  the  United  States.  In  addition,  a  multitude  of 
fixed  structures  ranging  in  size  from  small  aaymarks  to  offshore  towers 
are  used  to  "warn  the  navigator  of  dangers  and  oostruction. "  (U.S.  Coast 

Guard,  1972).  While  the  nature  of  these  fixed  aids  varies  widely,  ouoy 
classification  has  been  standardized  based  upon  their  maximum  di ampi-pr 
and  overall  length.  Standard  lighted  buoys  are  generally  made  of  steel 
plate  and  range  in  size  from  the  small  3  1/2  x  8  feet  (1.07  x  2.44  m.) 
buoy  to  the  large  9  x  32  feet  (2.74  x  9.75  m.)  shape.  Typically,  these 
buoys  are  held  on  station  by  a  rectangular  concrete  block  anchor  which  is 
connected  to  the  surface  floating  buoy  by  an  iron  cnain.  In  an  ideal 
situation,  the  length  of  this  mooring  chain  would  be  equivalent  to  the 
local  water  deptn.  Such  a  mooring  system,  however,  would  not  oe  very 
secure  as  any  horizontal  forces  resulting  from  either  wind,  waves  or 

current  would  quickly  move  the  buoy  off  station  or,  if  the  sinker  were 

large  enough,  cause  it  to  submerge.  Since  neither  of  these  alter¬ 

natives  is  desirable,  the  length  of  a  buoy's  mooring  chain  always  exceeds 
the  water  depth.  The  length  of  the  mooring  chain  is  dependent  upon  a 
host  of  factors  incluaing  the  average  water  depth,  the  maximum  current 
velocity  ana  the  size  of  the  buoy  itself.  The  ratio  between  the  mooring 
chain  length  and  the  water  depth,  or  scope,  has  been  found  to  generally 
be  about  2.5:1.  (Bitting,  1976).  There  are,  however,  certain  high 
current  areas  where  the  scope  necessary  to  securely  anchor  a  ouoy 
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is  considerably  larger  than  this  nominal  value.  Such  increases  in  scope 
result  in  a  subsequent  enlargement  of  a  buoy's  watch  circle  (the  area  on 
the  sea  surface  whose  perimeter  is  defined  oy  the  maximum  excursion  of  a 
buoy  from  its  anchor).  Aside  from  the  economic  considerations, 
increasing  a  scope  and  consequently  the  area  of  a  buoy's  watch  circle 
limits  the  effectiveness  of  the  aid  to  navigation  in  assisting  "the 
navigator  in  determining  his  position."  (U.S.  Coast  Guard,  1976). 

fixed  aids  are  preferred  to  buoys  because  of  the  positional  stabil¬ 
ity  they  provide  the  mariner.  Unfortunately ,  however,  floating  deoris 
and  ice  occasionally  strike  these  pilings  or  fixed  aids  with  sufficient 
force  to  damage  them.  Furthermore,  in  some  instances  the  lower  portion 
of  the  broken  pile  remains  submerged  while  still  implanted  in  the  seaoed 
thus  imposing  a  potential  hazard  to  vessel  traffic  in  the  area.  More 
frequently,  however,  these  pilings  are  struck  and  subsequently  broken  by 
barges  and  other  large  cargo  carrying  vessels.  That  such  accidents  take 
place  is  not  surprising  when  one  considers  the  rather  confined  quarters 
of  a  river  or  shipping  channel  as  compared  with  the  open  ocean,  and  that 
barges  are  frequently  rafted  together  around  one  tugboat  for  transport. 
The  effect  of  such  rafting  means  that  both  barges  and  tug  move  as  one 
massive  vessel.  Even  under  the  best  of  conditions,  the  maneuvering  of 
this  assemblage  is  difficult.  When  moving  downstream,  relative  to  the 
earth,  the  relative  flow  of  water  past  the  tug's  rudder  is  reduced  by 
the  velocity  of  tne  current  and  consequently  steerage  control  of  this 
tug-barge  system  becomes  all  the  more  difficult.  Considering  that  these 
piles  are  often  placed  in  locations  where  there  is  a  change  in  the 
direction  of  the  channel,  it  is  not  surprising  that  these  aids  are 
frequently  struck  or  grazed.  Besides  barge  traffic,  strong  currents 
have  a  tendency  to  scour  the  bottom  thereby  creating  the  possibility 
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that  it  will  erode  the  pile  where  it  enters  the  seabed. 

As  an  alternative,  an  articulated  spar  which  combines  the  most 
desirable  features  of  both  standard  buoys  and  conventional  fixea  aias 
might  be  used  as  a  functional  aid  to  navigation  in  these  situations. 
Such  an  aid  would  consist  of  an  anchor,  a  universal  joint  or  similar 
attachment  mechanism  and  a  rigid  member  extending  from  just  above  the 
channel  bottom  through  the  water  column  to  a  relatively  short  distance 
above  the  sea  surface.  Being  free  to  rotate  about  its  attachment  point, 
such  an  aid  would  tilt  when  grazed  by  a  passing  vessel.  In  addition, 
debris  or  flowing  ice  which  might  cause  a  pile  to  break  would  simply 
push  the  articulated  spar  aside  and  continue  to  move  downstream  or  else 
force  the  spar  underwater  while  passing  above  it.  Although  such  a  spar 
would  be  free  to  rotate  about  its  attachment  point,  its  rigidity  would 
restrict  its  watch  circle  to  a  very  small  area. 

1-2  Review  of  Relevant  Literature 

The  concept  of  an  articulated  spar  is  relatively  new  to  the  field  of 
Ocean  Engineering,  the  first  of  such  devices  being  the  product  of  oil 
company  research.  The  Elf-Ocean  Experimental  Oscillating  platform  is  a 
massive  cylindrical  articulating  structure  21  feet  (6. AO  m.j  in  diameter 
and  over  300  feet  (91.44  m.)  in  length.  (Villian,  1970}.  Originally 
placed  on  station  in  August  of  1968  in  the  Bay  of  Biscay,  this  platform 
was  designed  primarily  for  oil  exploration  purposes.  In  addition,  six 
large  cylindrical  subsurface  floats  provide  a  sufficient  amount  of  buoy¬ 
ancy  to  make  the  platform  stable  enough  for  general  offshore  research 
and  data  collection  purposes  as  well. 

Although  considerably  smaller  in  size,  another  articulated  tower 
which  had  its  origin  at  about  the  same  time  was  the  Single  Point  Mooring 
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(S.P.M.).  Also  the  product  of  oil  company  research,  these  moorings  were 
developed  for  the  loading  and  unloading  of  large  tankers  in  severe 
environments  ana  relatively  deep  water.  Such  moorings  were  designed  for 
the  purpose  of  being  a  much  less  expensive  and  more  easily  constructed 
alternative  to  the  deep  water  berthing  facilities  necessary  to  accommo¬ 
date  Very  Large  Crude  Carriers  (V.L.C.C.).  (Gruy  and  Kiely,  1977). 
While  somewhat  similar  in  concept  to  the  articulated  spar,  S.P.M.'s  also 
have  some  fundamental  differences.  The  first  S.P.M.  installed  at  Breqa, 
Libya  in  October  of  1969  for  the  Exxon  Company  serves  as  an  example. 
(Synodis  and  Flory,  1977).  This  S.P.M.  essentially  consists  of  an 
anchor,  a  universal  joint,  a  riser  section,  a  subsurface  buoyancy 
chamber,  a  mooring  float  and  some  connecting  chain.  In  addition,  one  or 
a  series  of  hose  sections  leads  from  the  swivel  joints  located  at  the 
top  of  the  buoyancy  chamber  to  either  the  vessel  being  loaded  or,  when 
not  in  use,  to  the  mooring  buoy  itself.  This  configuration  is  shown  in 
Figure  1-1.  From  this  sketch  it  can  be  seen  that  the  upper  portion  of 
the  S.P.M.  with  its  flexiole  anchor  chain  more  closely  resembles  a 
standard  mooring  system. 

In  the  early  1970's  a  series  of  light  beacons  known  as  Sarus  Towers 
were  constructed  and  installed  in  several  countries  includ¬ 

ing  Australia,  Papua,  New  Guinea  and  Indonesia.  The  Sarus  Tower  design 
combines  the  reliability  of  a  fixed  structure  with  the  low  capital  cost, 
ease  of  installation  and  flexibility  of  a  lighted  buoy.  Unlike  the 
articulated  spar,  the  main  body  of  the  Sarus  tower  is  composed  of  not 
one,  but  three  rigid  members  which  are  welded  together  end  co  end. 
While  the  exact  dimensions  of  these  towers  may  vary  from  site  to  site, 
the  installation  at  Port  Hedland,  on  the  western  coast  of  Australia 
serves  as  an  example.  (Whitaker,  1975).  At  this  location,  the  lower 
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Figure  1-1 

diagram  of  the  drega,  Libya  Single  Anchor  Leg  Mooring 


cylinder  is  a  42  incn  (106.68  cm.)  diameter  section  of  1/2  inch  (1.27 
cm.)  plate  steel  measuring  25  feet  (7.62  m.)  in  length.  The  upper 
cylinder  is  also  made  of  steel  plate  and  has  a  maximum  diameter  of  36 
inches  (91.44  cm.)  which  gradually  tapers  to  a  minimum  diameter  of 
approximately  24  inches  (60,96  cm.)  at  the  upper  end  of  its  40  feet 
(12.19  m.)  length.  A  large  subsurface  cylindrical  buoyancy  chamber  with 
dished  ends  serves  to  join  these  two  sections.  At  the  Port  Hedland 
installation,  the  buoyancy  chamber  has  a  diameter  of  9.5  feet  (2.90  m.) 
and  a  maximum  length  of  only  7.5  feet  (2.29  m.).  By  compartmentalizing 
this  buoyancy  chamber  into  four  sections,  one  section  can  become  flooded 
without  endangering  the  buoyancy  of  the  entire  structure. 

Another  group  of  articulated  light  beacons  are  those  located  at  the 
entrance  to  Genoa  harbor  in  Italy.  At  this  location,  the  nearly  one 
mile  long  entrance  channel  has  been  dredged  to  a  depth  of  approximately 
50  feet  (15.24  m. )  to  accommodate  large,  deep  draft  tankers. 
(Oell'Aggio,  1972).  However,  the  width  of  this  entrance  channel  ranges 
from  850  feet  (259.08  m.)  at  its  seaward  end  to  under  400  feet  (121.92 
m.)  at  the  jettyhead.  Consequently,  conventional  buoys  with  scopes  of 
about  2:1  were  found  to  be  unsatisfactory  for  defining  the  edges  of  this 
narrow  channel.  To  solve  this  problem  a  series  of  articulated  light 
beacons  were  designed  and  built  by  the  Resinex  Company  of  Iseo,  Italy. 
These  beacons  are  some  56  feet  (17.07  m.)  in  length  and  from  photographs 
appear  to  be  about  1  foot  (0.30  m.)  in  diameter.  They  are  composed  of 
an  anchor,  a  special  spherical  steel  joint,  a  galvanized  steel  pipe  and 
a  steel  frame  platform  which  is  mounted  at  the  surface  end  of  the 
structure  and  used  for  servicing  purposes.  Like  similar  structures,  a 
series  of  subsurface  floats  fixed  at  mid-depth  provide  additional  buoy¬ 
ancy  and  stability,  thus  reducing  oscillation  of  tne  structure  to  a  few 
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degrees.  Being  located  in  a  relatively  protected  environment,  these 
beacons  were  designed  td  withstand  the  combined  forces  of  approximately 
100  miles  per  hour  (160.93  km. /hr.)  winds  and  16  feet  (4.68  m.)  seas 
without  excessive  heeling. 

In  addition  to  being  used  as  stable  platforms  and  aids  to  naviga¬ 
tion,  an  articulated  spar  has  been  successfully  designed  and  constructed 
to  measure  oscillatory  flow  and  wave  direction.  (Lowe  et.  al.,  1974). 
Using  an  8.9  cm  diameter  air  filled,  filament  wound  pipe  attached 
through  a  universal  joint  to  a  bottom  anchor,  motion  was  measured  by  two 
orthogonally- mounted  accelerometers  located  in  the  top  portion  of  the 
spar.  Moored  in  33  feet  (10.06  m.)  of  water,  this  spar  behaved  as  a 
forced,  damped  oscillator.  Spectral  analysis  of  the  accelerometer 
records  showed  that  for  wave  periods  less  than  10  seconds,  the  tilting 
spar  showed  reasonable  agreement  with  a  computer  simulation  model  of  its 
behavior.  At  larger  wave  periods  the  model  predicted  tilt  angles  which 
were  too  low.  It  was  conjectured  that  this  was  the  result  of  the  com¬ 
puter  model's  use  of  constant  drag  and  inertia  coefficients. 

Excluding  Lowe's  tilting  spar  used  for  wave  measurements,  all  of  the 
aforementioned  devices  have  been  located  in  water  depths  of  40  feet  (12. 
19  m.)  or  greater.  To  the  author's  knowledge,  no  proposed  designs 
and/or  constructions  of  articulated  spars  in  the  15  to  30  feet  (4.57  to 
9.14  m.)  water  depths,  typical  of  river  channels,  currently  exist  in  the 
literature.  Without  exception,  the  greatest  design  force  imposed  upon 
any  of  these  existing  structures  is  the  result  of  wave  action.  In  the 
relatively  confined  area  of  a  river  or  shipping  channel  where  the 
proposed  articulated  spar  device  would  be  employed,  however,  current 
forces  are  expected  to  be  the  most  important  design  consideration. 


CHAPTER  II 


THE  ANALYTICAL  MODEL 

II-l  INTRODUCTION 

An  analytical  model  which  computes  the  total  force,  resultant  moment 
and  list  angle  of  a  spar  device  with  either  a  circular  or  elliptical 
cross-section  in  the  presence  of  wind,  current  and  wave  forces  has  been 
developed.  In  the  sections  that  follow,  a  detailed  derivation  of  this 
analytical  model  is  presented. 

I 1-2  Overview 

II-2.1  Theoretical  Considerations 

In  dealing  with  the  forces  on  a  three-dimensional  structure  such  as 
a  spar,  it  is  important  to  establish  a  coorainate  system  and  fix  its 
location  in  space.  The  origin  of  the  coordinate  system  used  in  this 
analysis  is  the  geometric  center  of  the  spar  cross-section  defined  by 
the  still  water  surface  at  zero  angle  of  inclination,  0,  as  shown  in 
Figure  II— 1.  The  velocity  components  U,V,  and  w  correspond  to 
velocities,  in  the  X,Y,  and  Z  directions  respectively  with  the  Z  axis 
being  defined  as  positive  upwards. 

A  combination  of  wind,  current  and  wave  forces  are  considered  to  be 
acting  on  the  articulated  spar  in  this  analysis.  Since  the  time 
responses  of  this  spar  system  to  changes  in  the  current  forces,  the  most 
critical  for  the  proposed  application,  are  shorter  than  the  time  rate  of 
change  of  the  current  field,  dynamic  effects  are  assumed  to  be  unimpor¬ 
tant  to  this  analysis.  Consequently,  this  analysis  assumes  the  articu¬ 
lated  spar  to  be  in  a  state  of  static  or  quasi-static  equilibrium. 

Fluid  particles  traveling  past  a  fixed  body  exert  both  normal  and 
tangential  forces  on  it.  Under  conditions  of  static  equilibrium,  the  * 


fh-  hydrodynamic  force 
fb-  buoyant  force 

Wgt  -  SPAR  WEIGHT 

0-  SPAR  ANGLE  OF  INCLINATION 

FIGURE  11-1  DEFINITION  SKETCH  OF  THE 
COORDINATE  SYSTEM 


body  must  remain  stationary,  requiring  that  the  summation  of  forces  and 
moments  must  be  zero: 


IF  =  0  (1) 

ZM  =  0  (2) 

where  o  may  be  any  point  on  the  body.  The  forces  on  the  spar  in  Figure 
II-l  may  be  described  in  vector  form  by 

IF  =  0  =  F^  ♦  Fd  +  wgt  (3) 

where  FH  is  the  hydrodynamic  force  composed  of  a  drag  and  inertia  com- 
ponent  following  the  approach  of  Morison,  et.  al.  (1950;,  F^  is  the 
buoyant  force  and  Wgt  is  the  spar  weight.  It  should  be  noted  that  the 
hydrodynamic  force,  also  includes  the  force  of  wind  on  the  exposed 
portion  of  the  spar.  Similarly,  the  moments  about  any  point  o  may  be 
described  by 

M0  ■  0  ■  «h  *  MB  *  V  (4) 

where  the  subscripts  are  the  same  as  those  used  in  Eq.  (3).  8y  applying 
these  two  fundamental  principles  of  statics  to  this  analysis,  the  list 
angle,  9,  of  the  articulated  spar  device  may  be  resolved.  A  further  term 
could  be  added  to  Eqs.  (3)  and  (a;  to  represent  a  lumped  mass  such  as 
might  occur  if  light  signal  equipment  ana  daymarks  were  added  to  the  spar 
or  to  denote  additional  floatation. 


In  this  analysis,  the  wind  velocity  is  assumed  to  act  on  that  portion 
of  the  spar  which  is  above  the  still  water  surface  to  induce  the  wind 
force  shown  in  Figure  II— 1 -  Similarly,  a  combination  of  current  and  wave 
particle  velocities  and  accelerations  are  assumed  to  act  on  the  submerged 
portion  of  the  spar  inducing  the  hydrodynamic  force  also  shown  in  Figure 
II— 1.  By  convention,  these  wind,  current  and  wave  particle  velocities 
and  accelerations  are  resolved  into  components  normal  and  tangential  to 
the  spar.  In  addition  to  these  two  normal  hydroovnamic  forces,  the 
buoyant  force  acts  vertically  upwards  from  the  center  of  buoyancy  and  the 
weight  acts  vertically  downward  from  the  center  of  gravity  as  described 
in  Equation  (3).  Once  these  forces  are  resolved,  moments  are  taken  about 
the  anchor  attachment  point.  In  the  convention  adoptea  in  this  analysis, 
the  overturning  moment  created  by  the  wind  and  water  particle  forces  acts 
clockwise  while  the  righting  moment  resulting  from  the  buoyant  force  acts 
counterclockwise  as  shown  in  Figure  11-1. 

The  total  hydrodynamic  force  on  the  spar  may  De  descriDed  by  tne 
Morison  equation  as 


fH  "  CD  *  0AcL  ^  S 


(5) 


where  the  first  term  on  the  right  hand  side  is  the  drag  force  and  the 
second  term  the  inertia  force.  In  the  drag  force  term  of  Equation  (5),  o 
is  the  fluid  density,  Ap  is  the  projected  area  normal  to  the  flow,  vn 
is  the  velocity  normal  to  the  object  and  Cg  is  a  dimensionless  drag 
coefficient.  The  absolute  signs  are  included  to  preserve  the  direction¬ 
ality  of  the  resulting  force  term.  The  inertia  portion  of  the  total 

hydrodynamic  force  is  composed  of  o,  the  density  of  the  fluid,  A  .. 

c  the 

,  „  dv 

cross-sectional  area  of  the  body,  l,  the  length  of  the  body,  Jt-  the 
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local  fluid  acceleration  and  CM  the  hydrodynamic  mass  coefficient. 


II-2.2  Steady  flow  forces 

Over  a  time  period  on  the  order  of  minutes,  those  forces  induced  by 
reasonably  steady  winds  and  currents  may  be  considered  steady  while  those 
induced  by  waves  will  be  time  dependent.  Consequently,  for  these  steady 

dv 

flow  forces,  the  ^ term  is  zero  and,  therefore,  the  total  wind  ana 
current  forces  may  be  completely  described  by  their  drag  component.  Thus 


Fwind  2  pAo^Vwind  lvwind 

n  n 


(6) 


current 


17) 


In  this  analysis,  one  may  specify  the  fluid  in  which  the  spar  is  immersed 
as  either  fresh  or  salt  water.  Although  the  density  of  brackish  river 
water  is  a  function  of  salinity  and  temperature,  the  variation  in  these 
properties  is  email  over  the  range  of  variation  expected  and  as  a  result, 
this  value  is  assumed  to  be  constant.  The  density  values  of  air  and 
water  selected  for  use  in  this  model  are  defined  at  a  temperature  of 
68°f  (20°C;. 

It  is  well  known  tnat  near  a  surface  such  as  land  or  water  the  wind 
velocity  profile  over  that  surface  is  not  uniform  but  varies  with  both 
time  and  altitude.  Ignoring  the  time  scale  variations,  the  wind  velocity 
profile  in  this  model  is  assumed  to  follow  the  well  known  Prandtl-von 
Karman  universal  velocity  distribution  law: 


U  In  (z/z  ) 

if  IH-TTo/TT 

o 


(8) 
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where  z  is  the  height  of  measurement  in  meters  and  zQ  is  an  effective 
roughness  parameter.  Since  only  large  velocity  winds  are  expected  to 
significantly  effect  a  spar's  list  angle,  a  zQ  value  of  .003  m.  as  sug¬ 
gested  by  Myers, et.al.  (1969)  for  strong  winds  is  used  in  this  model. 

As  a  result  of  this  assumed  distribution,  the  wind  speed  used  as 
input  to  this  computer  model  must  have  an  elevation  associated  with  it. 
8ased  upon  this  one  wina  data  point,  an  entire  wind  profile  including  the 
velocity  at  the  10  meter  level  is  determined.  Furthermore,  tnis  10  meter 
wind  velocity  value  is  used  as  input  to  the  Sverdrup,  Munk  and 
Bretschneider  (3.M.8.)  computations  described  below. 

11-2.3  Wave  Induced  Forces 

In  order  to  determine  the  total  water  particle  force  on  a  structure 
using  Equation  (5),  it  is  necessary  that  the  horizontal  and  vertical 
particle  velocities  and  accelerations  be  known.  Unlike  the  wind  and 
current  velocities  which  are  easily  measured,  the  wave  velocities  ana 
accelerations  required  for  this  analysis  are  often  estimated  using  simple 
wave  theory.  The  question  as  to  which  wave  theory  to  apply  to  determine 
the  appropriate  values  for  use  in  the  Morison  equation  has  been  studied 
extensively.  Patton  (1966)  states  that  for  a  body  that  is  small  compared 
to  the  wavelength,  the  water  wave  may  be  adequately  described  tv  Airy 
wave  thedry.  Dean  (1967)  has  shown  that  higher  order  Stokes  theories 
(3rd  and  5th)'  as  well  as  Cnoidal  wave  theory  are  not  uniformly  more 
apolicable  than  Airy  theory.  Further,  Dean  points  out  that  Airy  wave 
theory  satisfied  the  kinematic  boundary  conditions  at  the  surface  much 
better  than  these  higher  order  theories  for  shallow  water  waves  (Edge  and 
Myer,  1969).  In  general,  the  literature  seems  to  support  the  premise 
that  for  most  applications,  the  well  known  Airy  formulations  for  particle 


velocity  and  acceleration  are  the  most  logical  for  substitution  into  the 
Morison  equation.  Using  Airy  theory,  the  horizontal  and  vertical  compo¬ 
nents  of  the  local  fluid  particle  velocity,  if  and  it  respectively,  are: 


aak  cosh  k  (h  +  z) 

U  =~o  cosh  kh -  sin  <kx  "  at)  (?) 


= 


agk  sinh  k  (h  +  z) 
a  cosh  kh 


cos 


(kx  -  at) 


(10) 


where  g  is  the  gravitational  constant,  k  is  the  wave  number,  and  a  is  the 
wave  frequency.  The  local  fluid  acceleration  terms  in  the  horizontal  and 
vertical  directions  may  be  found  by  differentiating  Eqs.  (9)  and  (10) 
with  resoect  to  time,  t,  resulting  in 


3u 

7t 


-agk 


cosh  k  (h  -*•  z) 
cosh  kh 


cos  (kx  -  at) 


(11) 


3W 

TE 


=  -agk 


sinh  k  (h  +  z) 
cosh  kh 


sin  (kx  -at) 


(12) 


Equations  (9)  to  (12)  have  two  parts:  an  amplitude  term  and  3  phase 
term.  The  phase  term  in  each  describes  tne  oscillatory  nature  of  the 
wave  particle  velocity  and  acceleration.  The  amplitude  term  in  each 
equation  is  composed  of  a  constant  part  and  a  hyperbolic  part.  This 
hyperbolic  portion  accounts  for  the  exponential  decav  of  the  horizontal 
and  vertical  velocity  and  acceleration  components  with  increasing  water 
depth,  h  in  these  equations  is  the  time  dependent  total  depth  which  is 
defined  as 
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h  =  d  +  n  (13) 

where  d  is  the  still  water  depth  and  n  is  tne  corresponding  wave  height 
defined  by 


n 


2tt  t 

a  cos 


(14) 


where  a  is  the  wave  amplitude  and  T  the  wave  period. 

Because  of  the  nonlinearity  of  the  drag  coefficient  used  in  these 
drag  force  determinations,  the  total  water  particle  force  may  not  be 
determined  by  the  superposition  of  the  current  and  wave  forces  (Wu  and 
Tung,  1975).  Rather,  this  total  water  particle  force  must  oe  determined 
by  first  vectorally  summing  the  individual  velocity  components  and  then 
substituting  this  resultant  water  particle  velocity  into  the  drag  force 
term  of  the  Morison  equation.  Thus  tne  water  particle  orce  •’-ay  De 
defined  by 


t  1„,  |v  jv  Cn  +  oA  l  r> 
p—  =  2pAp  wpn  wpn  D  c  3E  C 


wo 


M 


(15.) 


where  the  water  particle  velocity,  v  ,  is  defined  as 


V  _  =  V  +  V 

wp  c  wa 


(16) 


where  v  is  the  current  velocity  and  v  is  the  water  particle 

(w  Wa 

velocity. 

In  this  analysis,  however,  no  consideration  is  given  to  the  possible 
interaction  effects  between  wave  induced  and  current  velocities.  Such 
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interaction  effects  have  been  examined  in  detail  by  Jonsson,  et.  al. 
(1970).  When  a  current  and  waves  are  traveling  in  the  same  direction, 
for  example,  wave  length  tends  to  increase  and  height  and  period  tend  to 
decrease.  These  effects  have  been  considered  significant  enough  that 
they  have  been  accounted  for  in  the  design  of  North  Sea  oil  drilling 
platforms  where  waves  are  the  dominant  forces  on  the  structures.  (Mes, 
1977).  However,  since  the  use  of  a  wave  correction  factor  to  account  for 
these  interaction  effects  necessitates  the  use  of  an  effective  average 
current  value,  this  formulation  was  not  included  in  the  analytical  model. 

The  waves  in  this  analysis  are  assumed  to  be  wind  generated  and 
consequently  act  in  the  same  direction  as  the  wind.  When  considering 
wave  forces,  input  data  is  eitner  supplied  in  terms  of  wave  height  and 
period  or  can  be  calculated  using  the  S.M.B.  method. 

II-2.4  Drag  Coefficient  Determination 

A  cylinder  with  a  circular  cross-section  is  the  most  common  struc¬ 
tural  element  in  the  ocean  and  as  a  result  a  great  many  experiments  botn 
in  the  laboratory  and  in  the  field  have  oeen  performed  in  an  effort  to 
define  precisely  the  drag  and  inertia  coefficients  of  the  Mo risen  equa¬ 
tion.  The  general  procedure  followed  in  these  experiments  has  been  to 
measure  the  forces  on  the  structure  as  well  as  the  wave  and  current  con¬ 
ditions.  Using  a  suitable  wave  theory,  the  water  particle  velocity  is 
subsequently  computed.  The  results  of  these  experiments,  however,  differ 
widely.  (Bretschneider  (1957),  Agerschou  and  Edens  (1965),  Wieselburger 
(1968),  Keulegan  and  Carpenter  (1956),  etc.).  Since  it  is  necessary  to 
employ  a  wave  theory  to  calculate  the  water  particle  kinematics  from  the 
measured  surface  profile,  the  resulting  range  in  these  coefficient  values 
reflects  not  only  possible  inaccuracies  in  the  Morison  formulation  but  in 
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the  wave  theory  used  as  well.  More  significantly,  it  was  discovered  that 
even  in  a  series  of  identical  waves,  scatter  was  found  in  the  wave  forces 
measured.  Further  analysis  has  shown  that  this  wiae  scatter  in  the 
published  data  is  related  not  only  to  the  wave  theory  selected  ana  the 
accuracy  of  the  measurement  system  but  also  to  the  roughness  of  the  pile, 
local  currents,  vibrations  in  the  test  piles,  the  proximity  of  neighbor¬ 
ing  piles  and  other  effects.  The  U.S.  Army  Coastal  Engineering  Research 
Center  has  examined  a  great  deal  of  this  data  in  an  attempt  to  simplify 
the  problem  of  selecting  the  proper  data  set  ana  have  subsequently  devel¬ 
oped  a  recommended  design  curve.  This  curve,  shown  in  Figure  IT.-2,  has 
been  incorporated  into  this  analytical  model  for  application  when  dealing 
with  circular  cylinaers.  One  should  note  for  purposes  of  comparison, 
Figure  II-2  also  presents  results  of  some  laboratory  and  field  measure¬ 
ments  of  drag  coefficients. 

In  Equation  (5),  it  was  shown  that  the  drag  coefficient  is  not  only  a 
function  of  Reynolds  number  but  roughness  and  length  effects  as  well. 

Further  it  was  pointed  out  that  data  regarding  these  effects  is  generally 
obtained  in  separate  analyzes  tb  establish  trends.  In  the  case  of  the 

dimensionless  spar  length-to-diameter  ratio  (L/D),  it  has  been  observed 
that  the  importance  of  this  parameter  increases  as  the  spar  length 
decreases.  In  such  cases,  the  pressure  at  the  end(s)  of  the  spar  is 
relieved  since  water  is  allowed  to  flow  around  the  end  of  the  spar  rather 
than  being  forced  to  one  side  of  it.  Consequently,  as  the  (L/D)  ratio 
decreases  so  does  the  drag  coefficient.  However,  empirical  data  regard¬ 
ing  this  effect  has  only  been  documented  for  Reynolds  numbers  in  the 

A  5 

range  10  to  10  .  Furthermore,  this  ratio  has  been  shown  to  be  most 
significant  for  small  [<  10]  (L/D)  ratios.  By  contrast  the  articulated 

spar  considered  in  the  analytical  model  is  generally  expected  to  have  an 
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(L/0)  ratio  larger  than  10  thus  minimizing  this  effect. 

In  general,  roughness  elements  tend  to  increase  the  effective  charac¬ 
teristic  dimension  of  a  body  and  thereby  effectively  increase  its 
Reynolds  number.  In  the  circular  cylinder  case,  decreased  roughness 
tends  to  result  in  the  transition  from  laminar  to  turbulent  flow  occurr¬ 
ing  at  lower  Reynolds  numbers.  In  addition  because  of  the  effectively 
increased  projected  area,  the  total  drag  force  on  such  a  cylinder  is 
larger  than  if  the  roughness  elements  were  absent.  While  such  effects 
have  been  measured  experimentally,  the  present  analysis  assumes  the  spar 
to  be  smooth  and  thus  ignores  these  effects. 

Drag  coefficient  information  on  an  elliptical  shape  is  by  contrast 
almost  non-existent.  Some  of  the  earliest  drag  coefficient  information 
on  this  shape  was  reported  by  Prandlt  (Daugnerity  and  Franzini,  1977). 
Unfortunately,  however,  this  information  only  examined  an  ellipse  whose 
major-to-minor  axis  ratio  was  4:1.  In  addition,  the  drag  coefficient  was 
determined  only  over  a  small  range  of  Reynolds  numbers  from  about  3  x 
104  to  3  x  105.  Hoemer  (1965),  after  examining  published  data  as 
well  as  performing  his  own  experiments,  developed  a  series  of  empirical 
relationships  for  the  drag  coefficient  on  an  ellipse  based  upon  Reynolds 
number  considerations  as  well  as  the  major-to-minor  axis  ratio  of  the 
ellipse.  For  flow  parallel  to  the  major  axis  these  relationships  were 
employed  in  the  analytical  model  to  compute  an  appropriate  drag 
coefficient. 

In  regard  to  flow  in  the  minor  axis  direction,  no  data  could  be  found 
in  the  literature  regarding  the  drag  coefficient  on  such  a  section.  How¬ 
ever,  it  was  noted  that  such  a  section  generally  resembles  a  flat  plate 
and  as  a  result  a  drag  coefficient  value  of  2.0  was  assigned  arbitrarily 
for  flows  in  this  direction.  While  such  a  value  may  be  rather  large,  in 


a  directionally  uniform  flow  field,  a  spar  would  align  itself  with  the 
current  and  thus  not  respond  to  these  off-axis  flows-  However,  if 
vertical  shear  were  present  in  the  water  column,  resolving  the  various 
velocities  acting  on  the  spar  into  components  results  in  only  a  small 
portion  of  any  off-axis  flow  being  actually  incorporated  into  the 
formulation  of  this  potentially  large  drag  coefficient  value. 
Furthermore,  since  any  off-axis  flow  would  also  generate  some  lift  on  the 
spar,  the  use  of  such  a  large  drag  coefficient  value  represents  a  poor 
attempt  at  including  this  effect.  Generally  the  lift  on  a  wing  section 
for  angles  of  attack  below  stall  is  proportional  to  the  sine  of  the  angle 
of  attack.  Unlike  a  wing  section, however,  the  elliptical  spar  section  is 
symmetric  about  its  minor  axis  and  has  a  blunt  trailing  edge  which  acts 
to  reduce  streamlining  effects.  Consequently,  the  coefficient  of  lift  of 
such  a  section  would  be  much  smaller  than  that  of  an  airfoil  shape. 
Therefore,  because  of  a  lack  of  adequate  lift  and  drag  information  as 
well  as  the  comparatively  small  likelihood  of  large  off-axis  flows,  a 
minor  axis  drag  coefficient  value  of  2.0  was  adopted. 

Nearly  all  of  the  previous  studies  which  define  the  drag  coefficient 
on  a  circular  cylinder  attempt  to  quantify  the  inertia  coefficient  as 
well.  Like  the  drag  coefficient  data,  a  great  deal  of  scatter  exists  in 
this  inertia  coefficient  information  also.  Unfortunately,  no  published 
data  or  empirical  relationships  could  be  found  relating  this  coefficient 
to  some  measurable  spar  parameter.  Consequently,  the  analytical  model 
employed  the  potential  flow  formulation  for  the  inertia  coefficient  of  an 
elliose: 


S*  *  (1  +  B/fl) 


(17) 
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where  A  and  B  are  the  major  and  minor  axis  lengths  of  the  ellipse, 
respectively.  Since  a  circle  represents  a  special  case  of  an  ellipse 
where  the  two  axes  are  of  equal  length,  this  formulation  reduces  to  the 
correct  potential  flow  value  of  2.0.  Generally  speaking,  potential  flow 
values  of  CM  tend  to  be  somewhat  higher  than  those  measured  in  the 
field.  Consequently,  the  forces  and  resulting  list  angles  computed  using 
these  potential  flow  values  will  tend  to  result  in  a  somewhat  conserva¬ 
tive  design  formulation. 

11-3  Computer  Model  Description 

The  forces,  moments  and  list  angle  of  an  articulated  spar  device  are 
calculated  according  to  the  procedure  described  above.  Using  Equations 
(6)  and  (15),  this  computer  model  consists  of  a  main  program,  ASB1,  and 
19  internally  called  subroutines.  Written  in  Fortran  IV,  this  program 
has  been  run  successfully  on  the  University  of  Rhode  Island's  ITEL  5 
digital  computer.  In  developing  this  model,  every  effort  has  been  made 
to  keep  this  program  as  general  as  possible,  and  as  a  result  the  program 
accepts  a  wide  variety  of  input.  The  large  number  of  arrays  which  may  be 
needed  by  this  program  require  that  320k  bytes  of  core  be  allocated  for 
its  execution. 

A  flowchart  summarizing  the  major  steps  in  the  program  is  shown  in 
Figure  II-3.  Initiation  of  this  computer  model  starts  by  defining  a 
series  of  fluid,  physical  and  environmental  parameters  whicn  are 
described  in  detail  in  the  users  manual  (Appendix  B).  In  summary,  the 
first  set  of  input  variables  serves  to  define  the  basic  fluid  properties 
of  viscosity  and  density  as  well  as  the  water  depth  in  which  the  spar  is 
located.  The  physical  dimensions  and  material  properties  of  the  spar 
including  the  orienting  vane,  if  present,  are  the  next  parameters  to  be 
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defined.  This  is  followed  by  a  description  of  the  distribution  of  the 
velocities  acting  on  the  spar.  Following  the  input  of  this  environmental 
information  a  series  of  subroutines  are  called  which  compute  information 
that  remains  unchanged  by  variations  in  the  spar’s  list  angle.  From 
these  subroutines,  sucn  information  as  the  weight  of  the  spar  in  air,  a 
variety  of  wave  information  and  the  wind  velocity  profile  are  obtained. 

Since  the  distribution  of  velocities  on  the  spar  may  have  large 
spatial  variations  in  the  vertical,  the  spar  is  divided  into  a  number  of 
elemental  areas.  8y  computing  the  total  force  of  each  of  these  areas  ana 
then  summing  over  the  entire  length,  the  total  force  acting  on  the  spar 
may  be  computed.  Since  the  variation  in  one  driving  force  profile  may  be 
much  larger  than  that  of  another,  three  sets  of  areas,  one  for  each 

force,  are  initially  defined.  However,  because  the  wave  and  current 
velocities  must  be  combined  in  order  to  determine  the  water  particle 
force,  this  number  of  area  sets  is  reduced  to  two;  one  for  the  wind 

forces  and  the  other  for  water  particle  forces.  In  this  reduction 

scheme,  the  maximum  number  of  equal  area  elements  describing  the  current 
and  wave  particle  velocities  is  selected. 

Having  defined  the  driving  force  velocity  profiles  acting  on  the 
spar,  its  orientation  relative  to  magnetic  north  is  ascertained  based 
upon  a  hierarchy  structure.  In  this  hierarchy,  the  force  due  to  current 
is  assumed  to  predominate  over  those  caused  by  wind  and  waves.  The 
coordinate  system  described  above  is  fixed  to  the  spar  and  as  a  result, 
the  bearing  of  the  buoy  relative  to  magnetic  north  can  oe  determined. 

The  computer  program  assumes  that  the  spar  will  orient  itself  such  that 
its  major  axis  is  parallel  to  the  direction  of  the  flow.  Further,  an 
orienting  vane  present  on  the  structure  will  always  act  as  the  trailing 
edge,  thereby  fixing  the  orientation  of  the  spar  relative  to  the  forces 


acting  on  it. 

Once  tne  orientation  and  number  of  eler.eni.al  areas  has  been  deter¬ 
mined,  the  relative  normal  and  tangential  velocities,  and  in  the  case  of 
waves,  accelerations  acting  on  each  of  these  areas  in  the  x  and  Y  direc¬ 
tions  are  determined.  Using  these  normal  velocities  to  first  compute  a 
series  of  elemental  Reynolds  numbers,  an  associated  drag  coefficient  for 
each  element  is  computed  using  one  of  the  previously  described  empirical 
relationships.  Subsequently,  the  force  on  each  of  these  elemental  areas 
is  determined  from  Equations  (6)  and  (15).  9y  summing  these  elemental 
forces,  the  total  wind  and  water  particle  force  acting  on  the  spar  is 
computed. 

These  driving  forces  produce  an  overturning  moment  whicn  is  computed 
next.  Buoyancy  forces  produce  a  righting  moment  which  acts  in  opposition 
to  this  overturning  moment.  Once  these  two  moments  have  been  calculated, 
they  are  compared  in  the  main  program.  Depending  upon  which  of  the  two 
moments  predominate,  the  list  angle  is  altered  and  the  computational 
process,  starting  with  the  determination  of  the  normal  and  tangential 
velocities  and  accelerations,  is  repeated.  When  the  two  moments  are 
equal  or  nearly  so,  depending  upon  the  convergence  criteria  specified, 
the  computation  is  complete.  At  this  time  the  wave  time  step  is  incre¬ 
mented  by  one  and  the  general  computational  process  is  repeated. 

When  the  forces,  moments  and  list  angles  have  been  determined  for  all 
of  the  wave  time  steps  then  one  final  subroutine  is  called.  This  subrou¬ 
tine  selects  an  appropriate  anchor  for  the  spar  based  upon  the  maximum 
net  buoyant  force  seen  by  the  spar  during  one  wave  period.  In  addition  a 
factor  of  safety  for  this  anchor  is  determined.  The  program  then  checks 
the  input  data  file  to  determine  if  there  are  additional  cases  to  be 
computed.  If  appropriate,  a  check  is  made  to  determine  if  during  the 


next  analysis  only  the  environmental  parameters  will  be  altered  or  if  the 
physical  dimensions  of  the  spar  will  also  be  changed.  By  returning  to 
the  appropriate  location,  the  new  data  is  input  to  the  main  program  ana 
the  computational  process  is  repeated.  When  no  further  data  is  supplied 
to  the  program,  it  terminates. 

11-4  Sensitivity  Analysis 

II-4.1  Introduction 

It  is  clear  that  the  analytical  model  which  has  been  developed  is  of 
a  very  general  nature  and  consequently  a  sensitivity  analysis  to  cover 
all  contingencies  could  not  be  made.  Therefore,  in  order  to  determine 
the  relative  importance  of  the  various  physical  and  environmental  param¬ 
eters  present  on  a  spar's  list  angle,  one  must  make  some  assumptions 
about  the  nature  of  these  possible  conditions.  The  first  step  in  the 
process  was  to  determine  the  significant  forces  involved  and  then  attempt 
to  ascertain  the  relative  importance  of  the  physical  and  environmental 
parameters  influencing  these  forces. 

II-4.2  Analysis  of  Forces 

The  analytical  model  previously  developed  allows  for  the  action  of 
wind,  wave  and  current  forces  to  act  upon  a  spar.  In  the  present  anal¬ 
ysis,  it  is  envisioned  that  an  articulated  spar  designed  using  this  model 
would  be  placed  in  the  relatively  confined  area  of  a  river  or  shipping 
channel  as  described  in  Chapter  I.  In  such  a  situation,  current  veloci¬ 
ties  could  be  on  the  order  of  2  to  3  knots  or  greater.  The  significance 
of  wind  forces  in  such  a  situation  is  best  illustrated  by  the  following 
example.  An  exceptionally  strong  wind  of  75  miles  per  hour  (120.70 
km. /hr.)  is  equivalent  to  a  velocity  of  110  ft/sec.  The  force  created  by 
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such  a  wind  velocity,  however,  is  smaller  than  that  created  by  a  current 
velocity  of  3  knots  (5.1  ft./sec.),  on  an  equal  area  basis.  For  example, 
on  a  one  square  foot  area,  with  an  assumed  drag  coefficient  of  1.0,  the 
wind  would  produce  a  force  of  28  pounds  (12.73  kg.)  while  that  created  oy 
the  3-knot  current  would  be  51  pounds  (23.18  kg.).  Aside  from  the  large 
velocity  differential,  the  difference  in  the  magnitude  of  these  two 
forces  is  the  result  of  density  differences:  that  of  water  being 
approximately  1000 times  greater  than  that  of  air.  More  significantly, 
while  such  current  velocities  are  expected  to  be  common  in  the  location 
where  such  a  spar  might  be  used,  the  wind  speed  used  above  is  representa¬ 
tive  of  hurricane  force  winds  and  thus  represents  an  extreme  case.  Fur¬ 
ther,  since  much  of  an  articulated  spar  is  assumed  to  be  submerged,  the 
area  on  which  the  current  force  acts  will  be  significantly  larger  than 
that  exDosed  to  wind  forces.  Consequently,  under  normal  conditions,  when 
small  daymarks  are  used,  it  is  reasonable  tb  assume  that  in  comparison 
with  current  forces,  wind  forces  may  be  neglected. 

In  general,  it  is  possible  to  omit  from  consideration,  the  action  of 
wave  forces  on  the  spar  as  well.  For  example,  consider  the  wave  forces 
created  by  a  large,  but  not  excessive,  wind  velocity  of  40  miles  per  hour 
(64.37  km. /hr.).  Such  a  velocity  will,  over  a  fetch  of  10,000  feet  (3.05 
Km.)  in  30  feet  (9.14  m. )  of  water,  produce  a  wave  with  a  significant 
height  of  2.5  feet  (0.76  m.)  and  significant  period  of  3.0  seconds  as 
predicted  by  the  S.M.3.  method.  Based  upon  Airy  theory,  such  a  wave 
would  result  in  a  maximum  horizontal  particle  ve’  city  of  2.6  feet  per 
second  (0.79  m./sec.).  While  this  value  is  approximately  50%  of  the  3 
knot  (1.55  m./sec.)  current  velocity,  several  factors  should  be  consid¬ 
ered.  First,  such  a  wind  velocity  and  fetch  condition,  while  not  rare, 
is  still  exceptional  in  the  semi-enclosed  areas  in  which  this  articulated 


spar  would  be  used.  More  importantly,  this  velocity  represents  a  maximum 
surface  condition  and  while  such  a  force  at  this  location  woula  produce  a 
sizable  moment,  this  force  is  exponentially  decreasing  with  both  time  and 
depth.  Consequently,  based  upon  the  scenario  developed  in  Chapter  I, 
combined  with  the  likelihood  of  forces  of  this  magnitude  occurring,  wave 
forces  on  a  spar  in  comparison  with  the  expected  current  forces  are 

generally  small  and  may  therefore  be  eliminated  from  further  considera¬ 
tion.  As  a  result  of  this  discussion,  current  forces  may  be  considered 
the  dominant  and  only  significant  force  acting  on  this  articulated  spar 
device. 

A  parametric  sensitivity  analysis  was  performed  on  many  of  the 
physical  and  environmental  parameters  which  effect  the  overall  list  angle 
of  the  spar.  As  this  analysis  proceeded,  it  became  apparent  that  because 
of  the  large  number  of  variables  involved,  different  spar  configurations 
might  result  in  very  different  sensitivity  results.  In  general,  by  vary¬ 
ing  one  parameter  while  keeping  all  others  constant,  trends  would  be 
established,  but  quantitative  assessment  of  the  various  differences  was 
limited  to  discussion  of  a  single  design.  Consequently,  the  parameters 
effecting  the  overall  list  angle  of  a  spar  buoy  are  most  meaningfully 
assessed  in  a  qualitative  rather  than  quantitative  sense.  Table  II-l  is 
a  list  of  a  variety  of  physical  and  environmental  parameters  with  a 
qualitative  assessment  of  their  influence  on  a  spar  buoy’s  overall  list 
angle.  These  influences  are  rated  as  either  significant,  moderate  or 
slight  and  are  based  upon  numerous  computer  simulations  under  a  multitude 
of  conditions.  Taole  II-l  serves  further  to  demonstrate  the  value  of 
this  analytical  model,  for  a  given  situation  where  many  of  these  param¬ 
eters  are  either  limited  or  fixed  due  to  certain  constraints,  the  user 
may  alter  the  remaining  parameters  wit:  an  whatever  physical  and 
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TABLE  II-l 


ANGLE 

Influence 
significant 

significant 
moderate 
moderate 
moderate 
slight 
slight 
slight 
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INFLUENCE  OF  VARIOUS  PARAMETERS  ON  A  SPAR'S  LiST 

Parameter 

current  velocity 

spar  weight  per  unit  length/buoyant 
weight  per  unit  length 

spar  length/ water  depth  ratio 

cross-sectional  area 

major  to  minor  axis  ratio 

wind  velocity 

wave  height 

anchor  attachment  offset 


28 


environmental  conditions  dictate  to  obtain  a  better  design  for  a  given 
situation.  Table  II— 1  may  therefore  be  used  as  a  guiae  as  to  the  effect 
of  varying  those  parameters  which  in  a  given  situation  may  be  altered. 

Of  the  parameters  presented  in  Table  II— 1 ,  the  current  velocity  and 
the  WPL/BPL  ratio  where  WPL  is  the  spar  weight  per  unit  length  and  BPL  is 
the  buoyant  force  per  unit  length,  were  judged  to  effect  significantly 
the  spar's  list  angle.  Figure  1 1-4  shows  the  effect  of  changes  in  the 
WPL/BPL  ratio  for  circular  cylinders  with  an  L/h  ratio  of  1.3,  where  L  is 
the  spar  length  and  h  is  the  local  water  depth.  At  low  velocities,  it 
can  be  seen  that  the  spar  is  very  sensitive  to  changes  in  the  WPL/BPL 
ratio,  while  at  large  velocities  this  effect  is  less  significant.  Total 
spar  submergence  occurs  at  a  list  angle  of  39.7°.  Visibility  consider¬ 
ations  dictate  an  upper  limit  on  WPL/BPL  of  approximately  0.65.  Spars 
with  ratios  greater  than  this  value  will  have  such  incrementally  small 
amounts  of  reserve  buoyancy  per  foot  that  their  list  angle  in  any  current 
would  be  excessive.  Material  properties  and  structural  considerations 
govern  the  lower  limit  of  this  ratio. 

Figure  II-5  shows  the  effect  of  various  L/h  ratios  on  the  list  angle 
of  a  circular  spar  in  various  uniform  current  fields.  It  can  be  seen  in 
this  figure  that  at  small  velocities  the  variation  in  L/h  is  most  signif¬ 
icant  while  at  large  velocities  a  spar's  list  angle  is  less  dependent  on 
this  ratio.  In  a  given  water  depth  as  the  spar  length  is  increased  so  is 
its  overturning  moment,  due  to  the  additional  weight.  To  compensate  for 
this  increased  overturning  moment,  the  spar  must  submerge  a  greater  por¬ 
tion  of  its  length.  Hence,  for  a  given  velocity,  spars  with  large  L/h 
ratios  exhibit  greater  list  angles  than  those  with  small  L/h  ratios. 
Therefore,  to  minimize  its  list  angle,  a  spar  shoulG  have  as  small  an  L/h 
ratio  as  possible.  However,  the  use  of  excessively  small  L/h  ratios 
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FIGURE  11-4  INFLUENCE  OF  CURRENT  ON  CIRCULAR 
CYLINDERS  OF  VARIOUS  WPL/BPL  RATIOS 

L/h  =  1.3 


LIST  ANGLE  (DEGREES) 


48 


L/h  =1.5 
L/h  =  1.4 


L/h  =1.3 


L/h  =1.2 


0.00  0.25 


0.50 


0.75 


1.00 


VELOCITY/ v/gD- 


1.25 


FIGURE  11-5  INFLUENCE  OF  CURRENT  VELOCITY  ON 
CIRCULAR  CYLINDERS  OF  VARIOUS  L/h  RATIOS 
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(i.e.,  approaching  1.0)  results  in  a  spar  becoming  submerged  sobner  than 
soars  with  large  values  of  L/h.  For  example,  in  Figure  11-5,  the  L/h  = 
1.2  line  terminates  at  a  dimensionless  velocity  value  of  0.75  and  the  L/h 
=  1.3  line  terminates  at  a  dimensionless  value  of  1.05,  indicating  the 
total  submergence  of  the  spar.  By  comparison,  spars  with  L/h  values  of 
1.4  and  1.5  become  totally  submerged  at  dimensionless  velocities  greater 
than  those  shown  in  this  figure. 

The  effect  of  changes  in  the  major-to-minor  axis  ratio,  A/B,  on  spars 
of  constant  projected  area  with  a  fixed  L/h  ratio  is  shown  in  Figure 
II-6.  It  can  be  seen  in  this  Figure  that  significant  reductions  in  list 
angle  can  be  achieved  by  cnanging  the  A/B  ration  from  1:1  (a  circular 
cross-section)  to  2:1.  Further  streamlining  of  the  elliptical  shape  from 
a  2:1  axis  ratio  to  a  3:1  ratio  results  in  less  significant  reductions  in 
spar  list  angle.  In  addition,  any  off-axis  flow  condition  would  be  sub¬ 
ject  to  large  projected  areas  at  these  large  A/B  ratios.  As  a  result  of 
this  analysis,  an  ellipse  with  a  2:1  A/B  ratio  was  selected  for  prototype 
testing. 
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FIGURE  11-6  INFLUENCE  OF  CURRENT  VELOCITY  ON 
VARIOUS  A/B  RATIOS  CONSTANT  PROJECTED  AREA 
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CHAPTER  III 


PROTOTYPE  TESTING 

III-l  Introduction 

The  validity  of  any  analytical  model  is  governed  by  its  ability  to 
describe  the  actual  situation  being  studied.  For  the  purpose  of  obtain¬ 
ing  validation  data  for  this  analytical  model,  a  series  of  laboratory 
tests  using  several  different  spar  shapes  were  performed  in  the  circulat¬ 
ing  water  channel  at  the  U.S.  Coast  Guard  Academy. 

III-2  Description  of  Facilities 

The  circulating  water  channel  used  in  these  tests  is  a  two-story 
high,  primarily  stainless  steel  structure  located  in  MacAllister  Hall  on 
the  grounds  of  the  U.S.  Coast  Guard  Academy  in  New  London,  Connecticut. 
Designed  at  the  University  of  Michigan  by  Vern  Phelps  and  Francis  ulgive, 
this  system  is  powered  by  an  axial  flow  propeller  pump.  By  adjusting  the 
propeller  speed,  velocities  up  to  a  maximum  of  10  ft./  sec.  (3.05 
m./sec.)  can  be  obtained  in  the  channel.  In  addition  to  being  able  to 
select  the  current  value,  plexiglas  viewing  ports  on  the  side  and  bottom 
of  the  circulating  water  channel's  10. 0  x  4.0  x  2.0  feet  (.305  x  1.22  x 
0.61  mj  of  usable  work  space  allow  the  user  to  observe  the  object  under 
study  from  many  different  angles.  In  addition,  colored  dye  is  availaole 
for  use  in  observing  general  flow  patterns  around  objects.  To  provide 
for  the  study  of  objects  which  must  be  anchored,  such  as  this  articulated 
spar,  a  series  of  screw  holes  are  located  in  several  places  on  the  bottom 
of  the  work  area,  witn  more  than  one  location  availade  at  a  given  time, 
for  comparison  purposes,  two  potential  spar  shapes  could  be  observeo 
together  in  the  same  flow.  Because  of  these  facilities  several  spar 
shapes  and  a  variety  of  orienting  devices  could  oe  testea  quickly  and 


easily. 

Current  velocity  measurements  in  the  circulating  water  channel  were 
made  using  a  Teledyne-Gurley  Moael  #625  6--cuo  pygmy  current  meter,  when 
in  use,  this  current  meter  is  positioned  on  a  steel  rod  which  is  also 
mounted  nn  a  cross-bar  across  the  tank  opening.  A  hand  held  revolution 
counter  connected  directly  to  this  steel  rod  allows  for  the  direct  deter¬ 
mination  of  the  average  flow  rate  over  any  known  time  interval. 

1II-3  Model  Descriptions 

Initially,  two  prototype  shapes,  one  with  an  elliptical  cross-section 
and  the  other  with  a  circular  cross-section,  were  constructed  for  testing 
purposes.  Later,  a  third  spar  with  a  more  streamlined  airfoil  shaped 
cross-section  was  built.  Figure  III-l  is  a  photograph  of  the  three  spar 
shapes  which  were  tested. 

Each  of  the  prototype  spars  was  constructed  from  a  styrofoam -type 
material.  Initially  a  template  of  the  desired  cross-section  was  con¬ 
structed  and  used  as  a  pattern  for  tracing  onto  the  foam  ooard.  Suose- 
quently,  thirty,  approximately  1  inch  (.54  cm.)  thick  foam  pieces  of  a 
given  cross-section  were  cut  and  then  fiberglassed  together  to  obtain  the 
desired  spar  shape.  Each  of  these  shapes  was  then  sandea  to  remove  any 
discontinuities,  and  a  bottom  attachment  plate  of  1/4  inch  (0.64  cm.) 
thick  plywood  was  glued  to  the  bottom  of  each  spar.  In  this  way,  ball 
joints  or  other  anchor  attachment  mechanisms  could  be  easily  secured  at 
any  desired  location  on  this  bottom  plate.  To  reduce  roughness  effects 
ana  minimize  drag,  each  model  was  then  painted  with  several  coats  of 
varnish. 

The  analytical  model  had  previously  shown  tnat  an  ellipse  with  a 
major-to-minor  axis  ratio  of  approximately  2 : 1  would  have  less  drag  than 
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Figure  III-l 


The  Three  Buoy  Shapes  Used  in  the  Circulating  Water  Channel 
Tes  ts 


a  circular  cylinder  of  the  same  cross-sectional  area  and  yet  be  less 
affected  by  off-axis  flows  than  ellipses  with  larger  axis  ratios.  Conse¬ 
quently,  the  first  model  which  was  constructed  had  a  major  axis  length  of 
4  inches  (10.16  cm.)  and  a  minor  axis  length  of  2  inches  *.5.08  cm.)  as 
shown  in  Figure  III-2A.  Such  dimensions  were  selected  to  approximate  a 
one-third  scale  model  of  the  prototype  articulated  spar.  In  order  that 
each  shape  have  the  same  buoyancy  per  unit  length,  the  subsequent  proto¬ 
type  shapes  were  constructed  so  that  their  cross-sectional  area  was  the 
same  as  that  of  the  elliptical  spar.  By  selecting  a  spar's  dimensions  in 
this  manner,  the  shape  with  the  smallest  list  angle  in  a  given  current 
would  be  the  one  with  the  least  drag  and  consequently  the  most  desirable 
for  field  test  purposes.  By  comparison,  the  circular  cross-sectionea 
spar  had  a  diameter  of  2.83  inches  (7.15  cm.)  as  shown  in  Figure  lii-20. 
After  some  initial  testing,  a  third  spar  with  a  cone  shaped  cross-section 
was  constructed  and  tested.  This  shape  had  a  chord  length  of  4.60  inches 
(11.68  cm.)  which  was  exactly  twice  its  maximum  width  of  2.30  inches 
(5.34  cm.).  This  spar  configuration  is  shown  in  Figure  III-2C.  To 
accommodate  the  2  feet  (0.61  m.)  water  depth  of  the  channel,  all  of  the 
spar  shapes  were  built  to  a  length  of  2  1/2  feet  (76.20  cm.).  Since  the 
full  scale  prototype  was  expected  to  be  considerably  greater  than  3  times 
this  length  (7  1/2  feet)  [2.29  m.  j,  geometric  similarity  between  these 
models  and  the  field  prototype  was  not  fully  maintained.  This  was  not 
considered  significant,  however,  because  the  models  were  constructed 
primarily  to  determine  an  optimized  spar  cross-section  and  as  such  their 
overall  length  was  of  secondary  importance. 

In  comparison  with  the  pressure  drag  forces  acting  on  the  spar,  how¬ 
ever,  the  wave-making  forces  created  by  the  spar's  piercing  the  surface 
were  minimal  and  consequently  all  velocity  scaling  considerations  were 
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,RE  III— 2  SCALE  DRAWINGS  OF  THE  THREE 
PROTOTYPE  SPAR  CROSS-SECTIONS 
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based  upon  the  Reynolds  number  rather  than  the  Froude  and/or  Weber 

numbers. 

III-4  Results 

The  initial  testing  was  done  using  only  the  elliptical  and  circular 
cross-sectioned  spars.  Initial  tests  with  the  circular  spar  showed  that 
at  any  significant  velocity  it  would  oscillate  in  a  direction  transverse 
to  the  current  flow.  This  type  of  motion  is  the  result  of  vortex  shed¬ 
ding  and  is  well  documented  in  the  literature  (Blevins,  1977;.  Tne 
period  of  this  oscillation  was  measured  and  generally  found  to  be  in  gooa 
agreement  with  that  predicted  by  empirical  formulations  based  on  the 
Strouhal  number.  The  Strouhal  number  is  generally  considered  to  be  a 
function  of  the  Reynold's  number  although  for  values  between  1.0  x  105 

and  3.0  x  105,  there  is  a  large  degree  of  spread  in  the  data  (Blevins, 
1978).  In  the  1.0  x  lO^1  to  8.8  x  10^  Reynold's  number  range  of  these 
tests,  however,  the  Strouhal  number  may  be  considered  constant  having  a 
value  of  0.19. 

Some  rather  unexpected  results  were  observed  when  the  elliptical 
cylinder  was  placed  in  the  circulating  water  channel.  It  was  anticipated 
that  tnis  elliptical  shape  when  anchored  from  a  point  near  its  leading 
edge  would  align  itself  such  that  its  major  axis  was  parallel  to  the 
direction  of  the  current  flow.  In  this  orientation,  an  ellipse  would 
oresent  a  very  small  projected  area  and  consequently  the  drag  force  on 
the  spar  would  be  minimized.  Unfortunately,  however,  using  this  ellip¬ 
tical  prototype  it  was  found  that  regardless  of  the  location  of  the 

ancnoring  point  on  the  bottom  face  of  the  spar,  this  shape  would  not 

align  itself  in  the  desired  direction.  Instead,  when  placed  in  a  current 
the  elliptical  spar  would  immediately  rotate  so  that  its  major  axis  was 


perpenoicular  to  the  direction  of  flow.  Once  aligned  in  this  position, 
the  ellipse  would  then  oscillate  transverse  to  the  flow  because  of  vortex 
snedding.  Milne-Thomson  (1938)  has  snown  that  an  ellipse  free  to  pivot 
about  its  geometric  center  will  rotate  until  the  two  stagnation  points  at 
either  end  of  the  major  axis  of  the  ellipse  lie  on  a  line  passing  through 
the  center  of  the  ellipse.  Therefore,  a  spar  mounted  in  this  manner 
would  be  expected  to  orient  itself  broadside  to  the  flow  direction.  How¬ 
ever,  this  elliptical  spar  was  attachea  to  the  oottom  very  near  its  leao- 
ing  edge,  thereby  greatly  reducing  the  magnitude  of  the  couple  which 
would  tend  to  rotate  the  buoy  model. 

Several  possible  causes  for  the  lack  of  directional  stability  in  tnis 
elliptical  spar  were  considered.  The  anchor  attachment  mechanism  used  in 
these  initial  tests  consisted  of  crudely  fabricated  pieces  of  aluminum 
which  in  principle  allowed  the  buoy  to  list  along  its  major  axis  and  to 
rotate  freely  aoout  its  vertical  axis.  Since  speculation  as  to  the 
reason  for  this  unusual  behavior  initially  centered  around  this  attacn- 
ment  device,  a  more  precise  mechanism;  a  small  ball  joint,  was  obtained 
and  used  in  all  future  testing.  Even  witn  this  device  secured  on  the 
bottom  face  of  the  spar  at  its  leading  edge,  the  elliptical  spar  turned 
broadside  to  the  current.  Having  eliminated  the  anchor  attachment  device 
from  suspicion,  other  possible  reasons  were  considered  to  explain  the 
spar's  unexpected  behavior.  It  was  hypothesized  or  deduced  that  compared 
with  an  airplane  wing  or  other  fixed  airfoil  shape,  the  additional 
degrees  of  freedom  of  this  system  were  responsible  for  its  lacK  of  direc¬ 
tional  stability.  Consequently,  either  a  perturbation  in  the  flow,  a 
discontinuity  in  the  actual  shape  of  the  soar,  a  slight  offset  in  the 
anchor  attachment  point  to  eitner  tne  left  or  right  of  the  centerline  or 
some  combination  of  these  possibilities  caused  the  initial  turning  of  tne 
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spar  in  one  direction.  Once  this  turning  began,  there  was  a  shift  in  the 
stagnation  points  and  a  differential  pressure  force  which  caused  a  net 
lift  force  to  begin  to  act  on  the  ellipse.  This  lift  force  would 
continue  to  increase  until  stall  was  reached.  However,  the  spar  woula 
continue  to  rotate  in  oppositon  to  the  increasing  drag  force  until  such 
time  as  the  stagnation  points  were  on  a  line  whicn  passed  through  the 
point  of  rotation.  This  stable  condition  would  thus  occur  when  the  spar 
was  broadside  to  the  current. 

Having  concluded  that  the  spar  alone  would  not  orient  itself  in  the 
desired  manner,  a  vane  made  of  sheet  aluminum  was  attachea  to  the  trail¬ 
ing  edge  of  the  spar  in  an  attempt  to  obtain  directional  stability. 
Because  of  the  problem  of  attaching  the  aluminum  to  the  foam  material  of 
the  buoy,  an  exceedingly  large  vane,  some  24  inches  (60.96  cm.)  in  width 
and  spanning  the  entire  length  or  the  spar  was  tried  initially.  In 
subsequent  tests  the  vane  was  shortened  by  cutting  from  the  trailing  edge 
what  was  deemed  to  be  excess.  The  results  of  this  series  of  experiments, 
also  performed  with  the  anchor  attachment  point  being  located  on  the 
bottom  face  of  the  spar  at  its  leading  edge,  showed  that  for  vane  widths 
greater  than  8  inches  (20.32  cm.)  (i.e.,  twice  the  buoy's  major  axis)  tne 
desired  effect  of  having  the  major  axis  of  the  ellipse  align  itself 
parallel  to  the  direction  of  flow  of  the  current  was  achieved.  Even  with 
these  large  vanes,  however,  the  spar  was  found  to  be  very  sensitive  to 
slight  perturbations  in  the  trailing  edge  of  the  vane.  Any  small  change 
in  this  trailing  edge  would  lead  to  the  entire  spar  assuming  a  list  angle 
which  was  transverse  to  the  flow.  However,  such  vane  widths  being  at 
least  twice  the  length  of  the  major  axis  of  the  spar  were  considered  to 
be  excessive  and  thus  this  design  was  excluded  from  further  consideration. 
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As  tne  vane  width  was  shortened  to  less  than  8  inches  (20.32  cm.;, 
the  overall  shape  of  the  spar  witr  a  blunt  leading  edge  and  a  sharp 
trailing  edge  began  to  resemble  that  of  a  wing  section.  The  consequence 
of  this  reduction  was  that  the  spar  became  increasingly  sensitive  to  very 
slight  perturbations  in  the  flow.  Once  this  perturbation  occurred, 
greater  lift  would  develop  on  one  side  of  the  spar  and  the  entire  spar, 
while  maintaining  its  orientation  to  the  flow,  would  assume  a  large 
transverse  list  angle  (Y  axis)  resulting  in  the  spar  being  completely 
suomerged.  While  it  was  at  times  possible  to  adjust  the  orienting  vane 

so  that  the  lift  forces  on  either  side  of  the  spar  were  nearly  equal  as 

shown  in  Fig.  Ill— 3 ,  any  slight  deflection  of  either  the  flow  in  front  of 
the  spar  cr  of  the  spar  itself  resulted  in  an  unstable  condition  from 
which  the  spar  never  recovered. 

Three  possible  mechanisms  to  reduce  this  lift  force  and  consequently 
improve  the  stability  of  the  spar  vane  system  were  tried  using  a  vane  4 

inches  (10.16  cm.)  in  width.  In  the  first,  a  hairy  fairing  was  addeo  to 

the  trailing  edge  of  the  vane.  With  this  device  no  notable  improvement 
in  system  stability  was  observed.  The  second  method  consisted  of  making 
a  series  of  horizontal  cuts  aoout  3/4  inches  (1.91  cm.)  deep  at  a  spacing 
of  approximately  1  inch  (2.54  cm.)  all  along  the  trailing  edge  of  the 
vane.  By  bending  these  resulting  tabs  in  alternating  directions,  a 
f lairing  of  the  trailing  edge  was  accomplished.  This  vane  modification 
improved  directional  stability  somewhat  but  not  without  greatly  increas¬ 
ing  the  drag  force  and  the  resulting  list  angle  as  well.  Furthermore,  it 
was  found  th?.c  only  slight  changes  in  these  tabs,  especially  those 
located  near  the  water  surface,  could  drastically  cnange  the  ouoy's 
stability.  Since  any  vessel  which  struck  such  a  spar  would  be  likely  to 
alter  these  trim  tabs,  especially  those  near  the  water  surface,  this  idea 


was  eliminated  frcm  further  consideration. 


Tne  Placement  of  a  series  of 


ncies  drilled  in  tre  vane  near  its  trailing  ecge  sic  appear  to  reduce  tne 
magnitude  of  tne  lift  force  somewnat ,  anc  consequently  increased  tne 
soar's  stability.  Unfortunately  nowever,  no  way  was  found  to  either 
significantly  reduce  or  stabilize  the  effect  of  this  lift  force. 

Measurements  of  the  list  angle  of  this  elliptical  spar  with  a  vane 
attached  were  taken  when  the  system  was  in  a  stable  mode.  These  measured 
values  were  found  to  be  in  good  agreement  with  tne  previously  develooed 
analytioal  model.  However,  the  analytical  model  (Chapter  Iii  die  not 
address  directional  stability  criteria  and,  therefore,  these  circulating 
water  channel  tests  were  essential  to  demonstrate  that  while  an  ellip¬ 
tical  spar  does  have  a  low  drag  shape,  directional  stability  Problems 
orecluce  its  use  without  a  more  in  dsptn  study  of  appropriate  orienting 
vanes . 

Since  the  ellipse  which  had  oeen  under  study  nao  a  major  to  minor 

axis  ratio  of  2:1  and  ellipses  of  larger  ratios  would  produce  even 

greater  amounts  of  lift  efforts  to  determine  a  staole  spar  focused  on  tne 

circular  cross-sectioned  soar.  The  oscillatory  motion  due  to  vortex 

shedding  nga  been  oreviously  observed  and  its  frequency  found  to  corre- 

S  v 

soona  with  that  predicted  by  the  Strouhal  number;  f  =  —  . 

It  is  well  known  that  one  of  the  most  common  ways  to  suppress  this 
shedding  motion  is  through  the  use  of  a  splitter  plate  or  vane  which  acts 
to  inhibit  the  formation  of  an  oscillating  pressure  field.  Consequently, 
~ne  of  these  elates  was  fashioned  from  a  s.neet  of  aluminum  sne  attached 
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returned  to  its  original  position. 

Having  met  with  this  success,  3  numoer  of  different  splitter  plate 
sizes  and  shapes  were  investigated  to  determine  the  minimum  size  vane 

which  would  ensure  the  stability  of  the  structure.  Tnis  minimum  size  was 

founa  to  be  about  1.25  diameters,  or  3.5  incnes  (8.89  cm.;  in  widtn,  and 

approximately  10.5  inches  (26.67  cm.;  in  length  (or  45 %  of  the  still 
water  submerged  lengtn  of  the  spar;.  Figure  Iii-4  is  a  photograpn  of 

this  design  being  tested  in  the  circulating  water  channel  at  a  velocity 
of  2.5  ft. /sec  (0.76  m./sec.;.  This  photograph  shows  the  vane  located 
near  the  center  of  the  overall  lengtn  of  the  spar.  Later  experimentation 
found  greater  directional  stanility,  as  determineG  by  the  rate  at  whicn 
the  soar  returned  to  an  eguilibrium  position  once  deflected,  could  oe 
acnieved  by  locating  the  vane  closer  to  the  top  of  the  spar.  while  tne 
ootimum  location  of  this  splitter  plate  relative  to  the  water  depth 
varies  with  the  current  velocity,  the  closer  the  top  of  the  vane  is  to 
the  water  surface  without  protruding  above  it,  the  better,  in  this  loca¬ 
tion,  the  splitter  plate  suppresses  the  oscillating  lift  force  wnic.n 
■would,  if  present,  produce  the  largest  overturning  moment  on  the  spar. 

Having  determined  the  minimum  size  splitter  plate  necessary  to  main¬ 
tain  the  stability  of  the  spar  structure,  numerous  measurements  of  tne 
current  velocity  and  the  corresponding  list  angle  were  mace.  This  data 
presented  in  Table  III— 1  was  tnen  compared  with  that  predicted  oy  tne 
analytical  model.  The  results  of  this  comparison  are  shown  in  Figure 
II 1—5 .  This  figure  snows  generally  good  agreement  between  the  observed 
and  predicted  results.  At  low  current  velocities,  the  theoretical  mocel 
aocears  to  overoredict  the  list  angle,  most  likely  because  of  tne  con¬ 
servative  nature  of  the  1.2  drag  coefficient  value  usee  in  tne  list  angle 
determination.  At  these  low  velocities  and  correspondingly  small  list 
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The  Circular  Spar  in  the  Circulating  Water  Channel  with  a 
Splitter  Plate 


lABLE  lil-l 

LIST  ANGlE  DATA  FOR  A  2.83  INCH  (,7.19  cm. )  DIhMEIER  SPAR,  3G  INCHES 
(76.20  cm;  IN  LENGTH  WITH  A  SPLITTER  PLATE 


Velocity 

FT. /SEC. ) 

[cm. /sec. ] 

MEASURED 

LIST  ANGLE 
(DEGREES) 

EST  .lMATEO 
ERROR 
(DEGREES) 

0.00 

[  0.00] 

-3.5 

+0.5° 

0.50 

[15.24] 

-1.0 

+0.  5° 

1.00 

[30.48] 

3.0 

+1.0° 

1.50 

[45.72] 

11.0 

+1.5° 

2.00 

[60.96] 

20.0 

+1.5° 

2.50 

[76.20] 

28.0 

+2.0° 

3.00 

[91.44] 

34.0 

+2.5° 

3.50 

[106.53] 

39.0 

+3.0° 

O  THEORETICAL  VALUES 


4.00 


{ 
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OBSERVED  VALUES  WITH 
ESTIMATED  ERRORS 


L 

-4.00  t— 

0.00 


0.80  1.60  2.40  3.20 

CURRENT  VELOCITY  (FEET/SEC) 


4.00 


FIGURE  111-5  CIRCULATING  WATER  CHANNEL.  THEORETICAL 
AND  OBSERVED  LIST  ANGLES  FOR  A  CIRCULAR  CYLINDER 


WPL/BPL  =  0.094 
L/h  =  1.25 
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angles,  the  length-to -diameter  ratio  of  the  soar  tnay  be  of  importance, 
thus  indicating  the  use  of  a  smaller  drag  coefficient  as  outlined  in 
Chapter  II.  At  higher  velocities  and  consequently  greater  list  angles, 
this  L/D  ratio  is  of  doubtful  importance.  However,  at  these  higher 
velocities,  the  surface  roughness  of  tne  spar  may  become  significant. 
The  effect  of  such  roughness  elements  is  to  increase  the  drag  coefficient 
and  consequently  tine  resulting  list  angle  would  be  greater  than  that 
predicted  by  theory. 

I I 1-5  Additional  Testing  Results 

In  an  effort  to  find  another  design  which  had  the  stability  charac¬ 
teristics  of  the  circular  cylinder  with  a  small  splitter  plate  attached 
yet  naving  a  smaller  drag  coefficient,  a  tniro  buoy  shape  was  constructed 

and  tested.  When  a  fluid  moves  past  a  cylinder  with  a  splitter  plate, 

two  identical  Foppl  vortices  form  on  either  sice  of  the  plate  in  the 

region  near  where  it  joins  the  cylinder.  The  location  of  these  vortices 

is  shown  in  Figure  II I -6.  Once  formed,  these  vortices  remain  stationary 

and,  therefore,  by  filling  in  these  areas  with  some  material  and  thus 
eliminating  them  altogether,  it  was  postulated  that  another  staole  buoy 
configuration  would  result.  The  resulting  buoy  had  a  conical  shaped 

cross-section  as  shown  in  Figure  III-2C.  While  this  streamlined  snape 
does  have  a  smaller  drag  coefficient  than  a  cylinder  of  tne  same  orojec- 
ted  area,  tests  performed  on  this  shape  showed  that  it  was  unstable  as  it 
sailed  back  and  forth  transverse  to  the  flow,  efforts  at  suppressing 
this  sailing  motion  met  with  only  limited  success.  A  hairy  fairing 

attacnea  to  the  trailing  edge  of  tnis  conical  snaoed  cevice  proved  to  oe 
tne  most  successful,  but  was  only  able  to  reduce  the  frequency  of  the 

sailing  motion  by  approximately  50%. 


mmsamm 
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WITH  A  SPLITTER  PLATE 


Lll-6  field  Model  Selection 


The  results  of  this  series  of  laboratory  tank  tests  demonstratea  tnat 
only  a  cylinder  with  a  circular  cross-section  and  a  splitter  plate  was 
directionally  staole  and  would  return  to  an  equilibrium  position  when 
temporarily  displaced.  In  audition,  generally  good  agreement  was  found 
between  the  measured  list  angle  of  this  spar  section  and  that  predicteo 
by  the  analytical  model.  Consequently  this  shape  was  chosen  for  the 
field  testing  phase  of  this  validation  study. 


CHAPTER  IV 
THE  F ItLO  TEST 


IV-1  Spar  Description 

A  full  scale  spar  prototype  was  constructed  for  the  purpose  of  vali¬ 
dating  the  analytical  model  under  actual  fielc  conditions.  ine  actual 
spar  was  constructed  from  18  feet  (5.49  m. )  of  6  inch  (15.95  cm.)  Giam- 
eter  schedule  40,  T6061  aluminum  pioe.  In  addition  to  this  pipe,  two  end 
caos,  a  mooring  attachment  device  and  a  vane  attachment  guide  were  faori- 
cated  3nc  then  welded  to  the  spar  section.  Drawings  used  in  the  con¬ 
struction  of  tnis  field  prototype  are  presented  in  Appendix  III.  While 
it  is  recognized  that  a  spar  of  these  dimensions  is  quite  small  in  ccm- 
oarison  with  standard  government  aids,  weighing  about  130  pounds  (59.1 
kg.),  this  spar  represented  the  most  practical  design  which  could  be 
easily  transported  and  deployed  without  mechanical  assistance.  Fabrica¬ 
tion  of  tnis  articulated  spar  was  performed  oy  the  Durant  Machine  Company 
of  Noank,  Connecticut. 

Once  faoricated,  a  vane  constructed  from  1/2  inch  (1.27  cm.)  plywood 
was  bolted  to  the  vane  attacnment  guide.  Tnis  vane  measured  11.72  feet 
(3.57  m.)  in  length  and  10.5  inches  (26.67  cm.)  in  widtn  and  was  situated 
so  that  the  oottom  edge  of  the  vane  was  A  feet  aoove  the  mooring  attacn¬ 
ment  point.  The  ratio  of  the  width  of  tnis  vane  to  the  diameter  of  the 
spar  was  about  1.58,  whicn  was  slightly  larger  than  the  minimum  1.25 
value  determined  in  the  circulating  water  channel  tests.  The  reason  for 
this  slight  increase  was  to  provide  seme  margin  of  safety  in  assuring 
that  the  soar  would  perform  as  exceeded. 

Once  the  vane  was  attached,  the  spar  was  painted  with  3  inen  ^ 20 . 32 
cm.)  alternating  orange  and  white  bands.  Figure  IV-1  is  a  photograph  of 
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the  completed  spar  device.  During  the  actual  deployment  a  guic*  flashing 
white  light  was  mounted  on  top  of  the  soar  as  an  additional  safer/ 

orecaution. 

A  swivel  joint  served  to  connect  the  spar  to  its  mooring.  Tnis 

swivel  could  be  attached  in  any  one  of  three  locations  along  the  ootrom 

of  the  spar  depending  upon  the  desired  angle  of  list  at  slack  current, 
this  attachment  swivel  and  the  three  possible  holes  are  shown  in  Figure 
IV-2. 

The  anchor  used  in  these  tests  consisted  of  four  100  pound  (45.^5 

kg.)  lead  blocks  placed  in  a  box  as  shown  in  Figure  IV- 3.  When  initially 
deployed,  each  of  these  four  blocks  was  individually  placed  in  the  cox, 
after  which  the  spar  was  attached  to  the  iron  eye  at  the  center  of  the 
box. 


IV-2  Siting  Considerations 

In  order  to  test  this  spar,  an  area  with  a  relatively  high  current 
and  shallow  water  near  the  University  of  Rhode  Island  was  requirea. 
Point  Cudith  Pond,  a  large  coastal  pond  located  along  the  southern  shore¬ 
line  of  Rhode  Island  was  selected  for  the  field  test.  it  nas 
oredominately  shallow  water  and  in  the  vicinity  of  its  narrow  opening  to 
the  ocean,  high  currents.  Figure  IV-4  is  a  map  of  the  pond  showing  the 
test  site  location.  While  currents  in  this  pond  are  strongest  in  the 
entrance,  safety  considerations  with  regard  to  the  large  amount  of  ooat 
traffic  here  precluded  placing  the  soar  in  this  area.  The  entrance 
channel  coming  into  the  pond  separates  into  two  sections  just  north  of 
the  narrow  entrance.  The  banks  of  the  much  less  traveled  right  hand  fork 
of  the  channel  is  occuoied  by  numerous  waterfront  cottages  and  the 
Department  of  Environmental  Management's  Marine  Experiment  Station.  At 


54 


Figure  IV-3 


The  Anchor  Used  for  Testing  the  Field  Prototype 


LOCATION 


i 

FIGURE  IV-4  MAP  OF  POINT  JUDITH  POND  AND  VICINITY 
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this  station,  two  piers  protrude  from  the  shore  to  the  edge  of  the 
channel,  making  an  excellent  observation  and  work  platform.  Further,  the 
availability  of  working  space  for  locating  the  instrumentation  recoraers, 
ample  shore  power  and  the  cooperation  and  assistance  of  the  employees  of 
the  Department  of  Environmental  Management  made  this  site  a  most  logical 
selection  from  which  to  perform  the  field  test.  The  widtn  of  the  channel 
at  this  location  is  very  narrow  being  only  about  200  ft.  (60.96  m. )  wide, 
and  consequently,  the  current  velocities  here,  while  not  as  strong  as 
those  at  the  pond's  opening  are  still  as  large  as  2.0  ft. /sec.  (60.96 
cm. /sec. ) . 

1V-3  Data  Collection 

The  spar  was  instrumented  with  two  Humphrey's  Model  CP17-0601-1 
inclinometers  mounted  orthogonally  to  one  another  on  the  inside  portion 
of  the  upoer  end  cap.  These  inclinometers  shown  with  the  end  cap  in 
Figure  IV -5  were  secured  to  the  spar  in  such  a  way  that  the  list  angles 
parallel  and  perpendicular  to  the  orienting  vane  and  consequently  the 
expected  flow,  were  measured.  A  multi-conductor  cade  passing  through  a 
water-tight  seal  in  the  end  cap  supplied  tne  necessary  input  voltage  and 
also  carried  the  output  signal  to  the  shore  based  recorders. 

The  inclinometers  used  in  this  field  test  consisted  of  a  series  of 
wire  wound  potentiometers  whose  resistance  was  controlled  by  the  motion 
of  a  pendulum  type  slider.  Thus  by  supplying  a  constant  current  at  a 
known  voltage  to  this  inclinometer  and  recording  the  output  voltage,  the 
angle  of  inclination  could  be  computed.  A  schematic  drawing  of  this 
inclinometer  system  is  shown  in  Figure  IV-6.  The  output  signal  from  ooth 
of  these  inclinometers  was  continuously  recorded  on  a  Hewlett-Packard 
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FIGURE  IV-6  INCLINOMETER  AND  LOAD  CELL  WIRING  DIAGRAM 


Model  710Q8  two-pen  strip  chart  recorder  using  different  colored  record¬ 
ing  pens  to  distinguish  the  two  cutout  signals.  The  input  signal  was 
periodically  monitored  to  ensure  its  consistency  with  a  digital  voltmeter. 

In  an  effort  to  measure  the  total  tangential  force  whicn  the  spar 
exerts  on  the  mooring  system,  a  load  cell  was  connected  between  the  moor¬ 
ing  eye  and  the  swivel  at  the  lower  end  of  the  spar.  3y  measuring  this 
force  and  oy  knowing  the  list  angles  in  ootn  the  X  and  Y  directions,  the 
total  force  acting  on  the  spar  could  oe  readily  computed.  Furthermore, 
oy  measuring  the  local  current  velocity  and  knowing  the  resultant  force, 
an  estimation  of  the  average  drag  coefficient  of  the  spar  could  be  mace. 
□Derating  in  much  the  same  manner  as  the  inclinometers,  changes  in  load 
were  decected  by  a  pressure  transmitter,  changes  in  loaG  were  detected  oy 
a  pressure  transmitter  which  converted  these  signals  to  changes  in  resis¬ 
tance.  9y  supplying  a  -xnown  voltage  to  this  load  cell  and  measuring  its 
cutout,  the  force  being  applied  could  be  computed  as  schematically  snown 
in  Figure  1V-6.  A  separate  cable  was  used  to  connect  this  piece  of 
instrumentation  with  the  necessary  shore  power  and  recording  oevice. 
Output  data  from  this  load  cell  was  recorded  on  a  single  channel  Linear 
Instruments  strip  chart  recorder.  Like  the  inclinometers,  the  output 
signal  from  this  load  cell  was  recorded  continuously  ana  the  input  signal 
was  periodically  cnecKed  with  tne  voltmeter  to  ensure  its  consistency. 

Current  velocity  measurements  in  the  vicinity  of  the  articulated  spar 
were  maue  using  a  Savonius  rotor- type  current  meter.  As  this  current 
meter  was  equipped  only  with  a  rotor,  only  speed  measurements  were  made 
Current  direction  was  recorded  based  upon  visual  observation  as  being 
either  ebo  or  flood.  These  current  soeeo  measurements,  taken  from  a 
small  boat  moored  immediately  adjacent  to  the  soar,  were  made  at  three 
different  locations  in  tne  water  column,  near  the  surface,  at  mib-oeptn. 
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and  at  the  bottom,  to  detect  any  vertical  variation  which  might  be 
present.  Because  of  the  rather  large  tidal  range  of  approximately  39 
inches  (99.06  cm.)  at  this  location,  no  attempt  was  made  to  determine  the 
exact  depths  at  which  these  current  measurements  were  taken.  Furtner, 
except  for  small  periods  of  time  near  slack  water,  little  or  no  vertical 
variation  was  found  in  these  current  profiles.  Because  of  the  time 
involved  in  making  these  current  measurements,  as  well  as  the  relative 
rate  of  change  of  the  current  speed,  these  readings  were  taken  at  approx¬ 
imately  30  minute  intervals.  Like  many  coastal  ponds,  a  large  amount  of 
seaweed  and  other  debris  was  present  in  the  water  coijmn.  Consequently 
the  problem  of  fouling  of  the  rotor  precluded  the  use  of  a  moored  current 
meter. 

The  water  depth  at  the  spar  was  another  parameter  which  was  measured 
as  part  of  this  field  test.  Using  a  stilling  well  mounted  on  the  Marine 
Experiment  Station  dock  nearest  the  soar,  relative  tidal  height  changes 
were  obtained.  By  measuring  the  depth  at  the  spar  and  simultaneously 
recording  the  stilling  well  reaaing,  a  correlation  between  other  stilling 
weil  observations  and  the  water  depth  at  the  spar  could  be  made.  Because 
of  the  large  tidal  range  at  this  location,  stilling  well  observations 
were  made  approximately  every  20  minutes  to  ensure  accuracy. 

IV-4  Results 

IV-4.1  Configuration  //I 

Oata  from  this  field  test  was  accumulated  over  the  period  21  May  to 
24  May  1979.  initially  the  soar  configuration  and  data  collection  proce¬ 
dure  described  above  were  used  and  those  results  are  presented  below  as 
Table  IV- 1.  The  time  recorded  in  this  table  has  as  its  origin  00:00  on 
the  21st  of  May.  Cursory  examination  of  this  data  indicates  that  the 
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TABLE  IV- 1 


CONFIGURATION  IP  1  FIELD  OATA 


MEASURED  * 


TIME 

VELOCITY 

DEPTH 

LIST  ANGlE 

(HOURS) 

(FT. /SEC.) 

(FEET) 

(DEGREES) 

16.50 

0.34 

11.68 

32.6 

17.00 

0.10 

11.91 

31.7 

17.50 

0.56 

11.90 

32.6 

18.00 

0.39 

11.77 

33.5 

16.62 

0.39 

11.48 

37.1 

23.50 

0.81 

8.39 

50 

24.20 

0.25 

8.35 

— 

24.83 

0.03 

8.50 

— 

25.50 

0.79 

8.78 

50 

26.12 

0.64 

9.46 

48 

L6.  75 

0.57 

9.75 

43 

27.33 

0.17 

10.25 

42 

27.75 

0.56 

10.57 

40 

28.42 

0.03 

11.04 

35 

29.00 

0.32 

11.33 

31 

29.52 

0.01 

11.48 

30 

30.00 

0.34 

11.42 

33 

30.67 

0.68 

11.16 

33 

31.17 

0.  8.1 

10.89 

35 

31.58 

1.00 

10.58 

35 

32.17 

1.52 

10.07 

42 

33.00 

1.69 

9.38 

47 

33  58 

1.69 

8.85 

50 

34.17 

1.40 

8.46 

50 

34.83 

1.42 

8.25 

52 

36.05 

0.32 

8.13 

50 

36.50 

0.  e>4 

8.21 

55 

*  Whole  number  values  from  hand  inclinometer 


magnitude  of  the  list  angle  of  the  spar  in  this  conf  a  .  ration  was  quire 
large  with  values  ranging  in  size  from  31.7°  to  55°.  Since  the 

inclinometers  mounted  in  the  spar  were  only  capable  of  resolving  angles 
up  to  45°.  much  of  the  inclination  data  presented  in  this  taole  was 

obtained  from  hand  inclinometer/protractor  measurement.  To  differentiate 
between  those  inclinometer  values  obtained  from  the  internal  inclin¬ 
ometers  and  those  measured  with  the  hand  hela  inclinometer/protractor, 
values  fitting  into  the  latter  category  are  listed  as  whole  numbers  in 

Table  IV- 1.  Figure  IV-7  is  a  photograph  of  the  spar  in  this  first 

configuration.  It  is  apparent  from  this  figure  that  a  great  deal  of  the 
soar's  18  feet  (5.49  m.)  of  length  is  above  the  water  surface  thereby 

explaining  the  large  list  angle  values.  In  additi.  n,  this  photograph 

shows  the  spar  oriented  with  the  splitter  plate  directed  into  the 

current.  Observations  of  this  unexpected  behavior  revealea  that  it  only 
occurred  at  very  low  velocities,  generally  less  than  0.25  ft. /sec.  \i.G2 
cm. /sec. ).  At  velocities  greater  than  this  threshold  value,  the  vane 

would  turn  so  that  it  was  no  longer  oriented  into  the  current  out  at  some 
angle.  As  the  velocity  continued  to  increase,  tne  spar  would  continue 
to  rotate  until  the  vane  was  situated  downstream  and  parallel  to  the  flew 
at  velocities  of  0.4  ft. /sec.  (12.19  cm. /sec.)  or  greater.  it  was 
initially  suspected  that  this  unusual  orientation  at  low  velocities  was 
the  result  of  the  buoyancy  of  the  vane  material.  However,  subsequent 
analysis  of  this  behavior  revealed  that  the  off-center  mooring  attachment 
point  was  the  cause.  In  this  first  configuration,  the  spar  was  anchored 
by  the  forward  most  of  the  three  anchor  attachment  holes  ki.e.,  the  hole 
furthest  from  the  vane). 

After  observing  the  behavior  of  the  soar  in  this  configuration  for 
approximately  two  tidal  cycles,  several  facts  became  evident.  Because  of 
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The  Field  Articulated  Spar  Anchored  in  Configuration  #1 


the  shallow  water  deptn  there  were  periods  of  time  ±n  which  the  water 
depth  was  less  tnan  half  of  che  total  spar  length.  Consequently,  wirn  so 
much  weight  in  the  air,  large  list  angles  were  encountered. 

This  experiment  was  extremely  useful  in  that  the  effects  of  spar 
submergence  compared  to  spar  length  became  readily  apparent.  It  was 
concluded  that  in  any  future  tests  the  total  spar  snould  be  shortened. 
Since  the  load  cell,  which  is  27  inches  (68.58  cm.)  in  length,  had  not 
provided  useful  data  during  this  test,  it  was  decided  to  remove  it  in 
future  field  tests. 

IV -4. 2  Configuration  42 

In  the  second  configuration,  the  load  cell  was  removed,  effectively 
decreasing  the  spar  length-to-water  depth  ratio,  and  the  mooring  attacn- 
ment  point  was  shifted  to  the  aftermost  of  the  three  attachment  holes 
(i.e.,  the  hole  nearest  the  vane)  to  prevent  any  reverse  orientation  of 
the  spar  at  low  velocities.  A  photograph  of  the  spar  in  this  new  config¬ 
uration  is  shown  in  Figure  IV -8.  Data  ootaineo  from  the  spar  in  this 
configuration  is  presented  in  Table  IV-2.  In  addition  to  the  list  angle 
values  measured  oy  the  inclinometer,  this  table  also  included  correspond¬ 
ing  outout  from  the  analytical  model.  Comparison  of  these  two  sets  of 
values,  however,  can  only  be  made  following  an  analysis  of  the  errors 
inherent  in  these  field  data  measurements  as  summarized  in  Table  IV-3. 

Slight  variations  in  the  actual  wind  and  current  profiles  as  well  3S 
vibrations  in  the  spar  contributed  to  the  uncertainty  in  the  inclinometer 
measurements  which  were  estimated  to  be  1.3°.  Current  and  water 
depth  measurements  were  taken  curing  this  field  test  for  use  as  input 
data  to  tne  analytical  model. 
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TA8LE  IV-2 


C0NFIGURAT1UN  #  2  FIELD  DATA 


VELOCITY 

DEPTH 

LIST  ANGLE 

CQMPUTtD 
LIST  ANGLE 

(FT./5EC. ) 

(FEET ) 

(DEGREES) 

(DEGREES) 

0.78 

10.46 

32.4 

34.2 

0.25 

11.43 

22.9 

25.8 

Q.  34 

12.37 

13.6 

2.0 

0.32 

12.17 

5.9 

1.6 

0.10 

13.40 

3.4 

1.0 

0.34 

13.17 

7.7 

1.6 

0.85 

12.60 

17.0 

15.2 

1.32 

12.25 

19.6 

24.2 

1.49 

11.94 

23.8 

27.2 

1.76 

10.40 

36.6 

38.1 

1.55 

9.83 

38.4 

39.8 

1.18 

9.62 

38.8 

40.4 

0.76 

9.58 

39.2 

40.2 

0.34 

9.59 

38.5 

39.9 

0.34 

9.76 

37.7 

38.9 

0.81 

9.94 

37.3 

38.6 

1.00 

10.34 

34.1 

35.8 

0.68 

,  10.78 

31.2 

32.3 

0.51 

11.20 

27.5 

27.3 

0.34 

11.74 

22.8 

20.9 

0.01 

12.20 

16  8 

1.3 

0.34 

12.62 

11.5 

1.8 

0.08 

12.82 

7.3 

1.1 

0.34 

12.85 

9.0 

1.7 

0.30 

12.37 

8.1 

1.7 

0.68 

12.57 

11.1 

14.3 

0.51 

12.34 

13.3 

14.3 

0.85 

12.09 

18.0 

20.9 

1.13 

11.77 

21.5 

26.4 

1.35 

11.32 

27.2 

31.2 

1.52 

10.88 

30.6 

33.7 

1.77 

10.21 

35.8 

38.9 

1.64 

9.77 

37.5 

40.2 

1.27 

9.54 

38.2 

41.0 

1.27 

9.54 

38.2 

41.0 

0.08 

9.47 

37.9 

40.5 

0.47 

9.67 

37.1 

39.5 

0.41 

9.91 

35.0 

38.5 

0.20 

10.57 

30.7 

32.9 

v 


68 


TABLE  IV- 3 


t 

FIELD  TEST  MEASUREMENT  ERRORS  i 


Variable 


Estimated 

Error 


List  Angle 


±1.3° 


Current  >10.76  cm. /sec. 

(Threshold  =  10.76  cm. /sec.) 


±10% 


Water  Depth 


±7.6  2  cm. 


i 
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A  3a von i us  rotor,  which  is  well  known  in  tne  literature  for  its 
non-linearity,  was  used  to  determine  current  profiles  in  the  vicinity  of 
the  spar.  Typically  these  current  meters  have  a  non-zero  threshold,  a 
tendency  to  stall  at  very  high  speeds,  and  exhibit  more  rapid  accelera¬ 
tion  than  deceleration  (Seardsley ,  et.  al.,  1S67).  In  addition  tilt 
effects  have  been  shown  to  significantly  effect  the  accuracy  of  measure¬ 
ments  made  with  these  devices  although  their  exact  impact  varies  from  one 
instrument  to  another  as  well  as  over  various  speed  ranges.  During  the 
field  test,  however,  it  was  observed  that  at  low  velocities  (less  than 
0.35  ft. /sec.  (10.67  cm. /sec.))  the  current  meter  readings  were  subject 
to  large  fluctuations  resulting  from  the  inability  of  the  magnetic  reed 
switches  in  the  current  meter  to  function  properly.  Consequently ,  tne 
errors  associated  with  these  low  velocity  measurements  were  extremely 
large  and,  therefore,  such  values  were  omitted  from  furtner  study. 
Current  values  greater  than  0.35  ft. /sec.  (10.67  cm. /sec.)  were  estimatec 
to  be  accurate  to  within  _+  10%. 

The  water  depth  measur  .ent  used  as  an  input  variable  to  tne  analyt¬ 
ical  model  is  the  distance  from  mean  sea  level  to  the  oottom  of  tne 
spar.  These  water  depth  values  were  estimated  from  stilling  well 
observations  and  a  single  correlation  measurement  of  the  total  water 
depth.  Although  the  stilling  well  observations  were  virtually  error- 
free,  the  total  depth  determination  was  subject  to  an  error  estimatec  to 
be  +  2  inches  (+_  5.08  cm.).  In  addition,  the  distance  from  the  oase  of 
the  soar  to  the  pond  bottom  was  measured  to  3n  accuracy  of  +  1  inch  (+_ 
2.5u  cm.).  Consequently,  the  total  error  associated  with  these  water 
deoth  measurements  was  +_  3  inches  (_+  7.62  cm.). 

A  sensitivity  analysis  was  performed  to  determine  the  significance  of 
these  current  velocity  and  water  depth  uncertainties  on  predicting  th~ 
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list  angle  of  the  spar.  Assuming  a  +_  10%  error  in  current  values  greater 
than  0.35  ft. /sec.  (10.67  cm. /sec.),  the  change  in  the  spar's  list  angle 
was  found  to  be  on  the  order  of  0.1°  to  0.2°.  Consequently,  for  this 
buoy  configuration,  potentially  large  errors  in  the  current  velocity 
measurements  did  not  significantly  effect  the  spar's  list  angle. 

Figure  IV-9  shows  the  results  of  a  sensitivity  analysis  of  the  effect 
of  changes  in  the  water  depth  over  the  range  experienced  in  this  field 
test  on  tne  spar's  list  angle  assuming  three  separate  current  veloci¬ 
ties.  The  relative  steepness  of  the  slope  of  these  lines  indicates  that 
this  field  test  soar  buoy  was  very  sensitive  to  slight  changes  in  the 
water  depth  especially  at  low  velocity  values.  Therefore,  of  the  two 
parameters  measured  for  use  as  input  data  to  the  analytical  model,  only 
errors  in  tne  water  depth  were  found  to  significantly  affect  tne  spar's 
predicted  list  angle. 

Figure  IV-10  is  a  parity  plot  of  the  measured  vs.  computed  list  angle 
data  with  error  bars  indicating  the  uncertainties  associated  with  indi¬ 
vidual  water  depth  and  inclinometer  determinations.  The  relative  varia¬ 
tion  in  the  width  of  these  error  bars  results  from  absolute  variations  in 
the  current  velocity  and  water  depth  measurements  from  which  these  list 
angle  values  were  computed.  Comparison  of  model  predictions  ana  field 
measurements  as  noted  in  Figure  IV-10  snows  excellent  agreement  with  a 
correlation  coefficient  of  0.986. 

IV-A.3  Splitter  Plate  Removal 

As  the  final  phase  of  this  field  test  the  splitter  plate  was  removed 
from  the  soar  to  determine  the  effect  such  a  change  would  have  on  the 
performance  of  the  spar.  It  was  expected  that  in  stronger  current  situa¬ 
tions  (>0.5  knots)  [>25.76  cm. /sec.]  that  some  vortex  sheading  motion 
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would  be  observed.  Equipment  problems  prevented  the  accurate 
determination  of  the  spar  list  angle,  but  no  unusual  list  angles  were 
observed  when  the  splitter  plate  was  removed.  However,  the  expected 
vortex  snadding  motion  did  not  occur. 

As  was  true  in  configuration  #2,  the  spar  was  moored  at  the  aftermost 
of  the  three  mooring  attachment  points.  It  was  observed  that  turning  the 
spar  about  its  vertical  axis  did  not  result  in  an  arbitrary  reorientation 
of  the  spar.  Rather,  the  spar  always  tended  to  return  to  a  position 
where  the  mooring  attachment  plate  was  parallel  to  the  current. 
Initially,  it  was  speculated  that  the  welded,  vane  attachment  point  was 
acting  as  a  splitter  plate  and  thus  preventing  any  vortex  shedding 
motion.  The  hypothesis  was  rejected  because  subsequent  analysis  revealed 
that  the  uncentered  mooring  attachment  resulted  in  a  torque  being  placed 
on  the  spar.  This  torque  combined  with  the  asymmetry  of  the  spar  result¬ 
ing  from  the  vane  attachment  plate  resulted  not  only  in  the  spar  orient¬ 
ing  itself  into  the  current,  but  prevented  vortex  shedding  from  occurring 
at  these  velocities  by  inhibiting  transverse  motion. 


CHAPTER  V 


SUMMARY 

V-l  Review  of  Results 

An  articulated  spar  as  descriDed  in  this  study  mignt  serve  as  a  func¬ 
tional  alternative  to  buoys  ana  fixed  aids  to  navigation  which  are 
presently  used  in  rivers  and  other  high  current  areas.  When  struck, 
fixed  structures  are  frequently  destroyed,  whereas  an  articulated  spar, 
being  free  to  rotate  about  its  attachment  point,  would  either  move  out  of 
the  way  or  be  forced  underwater  until  the  object  passed  over. 

To  investigate  such  a  structure,  an  analytical  model  capable  of 
resolving  the  forces,  moments  and  list  angle  of  an  articulated  spar  oeing 
acted  upon  by  the  forces  of  wind,  current  and  waves  was  developed.  An 
analysis  of  the  major  physical  and  environmental  parameters  influencing 
these  forces  and  moments  was  made  subsequently.  In  the  relatively  con¬ 
fined  areas  of  a  river  or  shipping  channel,  only  current  forces  were 
found  to  contribute  significantly  to  the  motion  of  the  spar.  Further 
analysis  revealed  that  in  addition  to  current,  only  one  other  parameter, 
the  spar  weight-per-unit-length/buoyant-force-per-unit-length  ratio, 
affected  significantly  a  spar's  list  angle.  An  investigation  of  several 
shape  related  parameters,  including  the  spar  length -to-diameter  ratio  ana 
the  effect  of  varying  cross-sectional  areas,  demonstrated  that  these 
variables  affect  moderately  the  determination  of  a  soar's  list  angle. 

To  validate  this  analytical  model  a  series  of  tests  were  performed  in 
a  circulating  water  channel.  An  elliptically  cross-sectioned  spar  with  a 
major-to-minor  axis  ratio  of  2:1,  moored  near  its  leading  edge,  was  found 
to  be  unstable  in  a  current  flow.  Since  the  analytical  model  did  not 
address  directional  stability,  this  result  was  not  predicted.  Rather 
than  align  itself  with  its  major  axis  parallel  to  the  current,  the  spar 
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would  rotate  so  that  this  axis  was  perpendicular  to  the  flow.  Careful 
study  revealed  that  this  oenavior  was  the  result  of  an  asymmetric  lift 
force  being  generated  either  by  slight  perturbations  in  the  flow  or 
slight  asymmetries  in  the  spar  shape.  Subsequent  attempts  at  improving 
the  directional  stability  of  this  elliptical  spar  met  with  only  very 
limited  success.  As  a  result,  the  scope  of  the  testing  program  was 
expanded  to  include  identifying  a  spar  cross-section  which  was  direction- 
ally  stable.  Numerous  tests  of  both  circular  and  tear-drop  cross- 
sectioned  spars  revealed  that  of  those  tested  only  a  circular  cross- 
sectioned  spar  with  a  splitter  plate  had  the  desired  directional  stabil¬ 
ity  characteristics  over  a  broad  range  of  flow  conditions.  A  study  of 
different  splitter  plates  found  that  the  minimum  size  necessary  to  ensure 
stability  was  1.25  spar  diameters  in  width  and  approximately  45&  of  tne 
water  depth  in  length.  Comparison  of  the  performance  of  this  design 
shape  with  that  predicted  by  the  analytical  model  showed  generally  good 
ag reement. 

Based  upon  the  results  obtained  in  the  circulating  water  cnannel 
tests,  a  field  prototype  spar  was  faoricated  from  nominal  6  inch  (.16.33 
cm.;  diameter  schedule  40  aluminum  pipe.  Like  the  circulating  water 
channel  prototype,  this  spar  was  equipped  with  a  splitter  plate  as  well 
as  two  orthogonally- mounted  inclinometers  to  measure  list  angles.  Test¬ 
ing  of  this  prototype  spar  took  place  over  a  four-day  period  during  which 
current  velocity,  water  depth,  list  angle  and  axial  load  data  were 
recorded.  After  one  day  of  testing  tne  initial  buoy  configuration  was 
altered.  Shallower  than  expected  water  depths  at  the  test  site  resulted 
in  spar  list  angles  which  were  frequently  larger  than  the  45°  limit  of 
range  of  the  inclinometers.  Since  it  was  not  possible  to  mcve  the  spar, 
the  original  configuration  was  modified.  By  removing  the  load  cell,  from 
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which  no  meaningful  data  had  been  obtained,  the  water  depth  at  the  spar 
was  effectively  increased  by  27  incnes  (0.61  m.). 

In  the  first  configuration,  it  was  observed  that  the  spar  would 
orient  itself  so  that  the  splitter  plate  was  upstream  of  the  spar  at  low 
current  velocities  (  <0.25  ft. /sec.)  [7.62  cm. /sec.].  To  prevent  this 
undesirable  orientation  from  reoccurring,  the  mooring  attachment  point 
was  shifted  from  the  forwardmost  to  the  aftermost  location  (i.e.,  to  the 
hole  nearest  the  splitter  plate).  Before  comparing  the  results  obtained 
from  the  spar  in  this  second  configuration  with  those  predicted  by  the 
analytical  model  an  analysis  of  the  errors  inherent  in  these  field 
measurements  was  made.  This  analysis  revealea  that  the  error  associated 
with  tne  inclinometers  was  on  the  order  of  _+  1.3°.  Current  velocity 
and  water  depth  measurements  taken  during  this  field  test  were  used  as 
inDut  data  to  the  analytical  model.  An  estimated  error  of  +  10%  in  the 
measured  current  velocity  was  found  to  alter  a  spar's  predicted  list 
angle  by  only  0.1°  to  0.2°.  Depth  measurements  were  estimated  to 
h3ve  an  accuracy  of  +  3  inches  (  _+  7.62  cm.).  Over  the  range  of  water 
depths  experienced  in  this  field  test,  uncertainties  of  this  magnitude 
were  found  to  significantly  affect  the  list  angle  predicted  by  the 
analytical  model.  Oespite  these  inherent  errors,  comparison  of  tne 
analytical  model  predictions  with  the  field  test  measurements  snoweo 
excellent  agreement,  having  a  correlation  coefficient  of  0.986. 

V-2  Conclusions 

The  general  conclusion  of  this  study  is  that  the  development  of  an 
analytical  model  of  the  forces,  moments  and  resulting  list  angle  on  an 
articulated  spar  of  circular  cross-section  was  successful.  Validation  of 


this  model  for  the  test  spar  over  a  broad  range  of  environmental  condi¬ 
tions  showed  generally  gooo  agreement  between  preaicted  and  measured  list 
angles.  In  the  field  test  where  there  were  several  sources  of  error,  the 
results  obtained  showed  excellent  agreement.  Therefore,  this  analytical 
model  would  serve  as  one  component  of  a  design  tool  to  predict  the  per¬ 
formance  of  a  variety  of  articulated  spars  of  circular  cross-section 
under  a  range  of  current  velocities.  Tne  use  of  this  analytical  model 
for  predicting  the  list  angle  of  elliptically  cross-sectioned  spars  is 
not  recommended  without  aaditional  consideration  of  directional  stability. 

V-3  Suggestions  for  Future  Study 

During  tne  course  of  this  work,  it  became  evident  that  research  in 
several  areas  could  yield  additional  information  which  woulo  improve  the 
design  of  the  articulated  spar. 

V-3.1  A  Better  Shape 

The  use  of  the  circulating  water  channel  in  tnis  experiment  was 
invaluahle  in  facilitating  the  testing  of  numerous  spar/vane  shapes  under 
a  variety  of  conditions.  Of  the  sections  examined,  the  circular  shape 
used  as  the  basis  for  this  study,  was  unquestionably  the  most  staple. 
However,  for  a  portion  of  this  project  which  was  not  scheduled,  the  test¬ 
ing  with  these  facilities  was  constrained  by  time  and  consequently  was  by 
no  means  exnaustive.  Therefore,  one  potential  area  of  future  research 
would  be  the  examination  of  additional  low-arag  but  airectionally  stable 
shapes.  Despite  tne  problems  with  spar  length,  it  is  felt  that  because 
the  ccmputer  analysis  of  this  spar  incorporates  the  use  of  elemental 
areas,  Reynolds  scaling  of  any  shape  could,  with  some  slight  modifica¬ 
tions  to  the  drag  coefficient  subroutines,  be  accomplished  for  any  spar 


shape.  Unless  such  a  shape  is  commercially  available,  the  field  testing 
of  sucn  a  design  could  require  a  special  extrusion  ana  thus  prove  to  be 
costly. 

One  possiole  suggestion  for  such  a  shape  would  be  a  hybrid  system 

involving  a  wing- shaped  device  whicn  is  clampeu  over  a  circular  spar.  By 
oeing  free  to  rotate  about  the  spar  axis,  these  airfoil  shapes  should 

reduce  the  overall  drag  on  the  system  and  be  free  to  align  themselves 

with  the  prevailing  currents.  Presently,  similar-in-design  but  much  more 
massive  clamshell  structures  are  oeing  successfully  usea  on  oil  drill 
risers  to  reduce  the  drag  due  to  current  forces  (Chaaakoff,  1978;. 

V-3.2  field  Test  Modifications 

It  was  clear  from  the  field  test  that  for  the  conditions  and  equip¬ 
ment  available  some  very  useful  information  was  gathered.  however, 
several  important  improvements  could  be  made  in  future  field  tests. 

Since  the  current  velocity  profile  and  water  depth  are  critical  inout 
parameters  to  the  analytical  model,  a  more  accurate  determination  of 
these  values  would  lead  to  more  precise  predictions  of  the  list  angle. 

Unfortunately  most  areas  with  strong  currents  generally  tend  to  nave 
large  amounts  of  bulky  material  sucn  as  seaweed  suspended  in  the  water 
column.  Consequently,  a  permanently  mounted  current  meter  is  likely  to 
soon  become  fouled  and  thus  its  use  is  not  suggested.  However,  a  perma¬ 
nent  anchor  with  a  wire  attached  might  be  deployed  in  such  a  way  that  a 
current  meter  could  be  easily  clampea  to  it  and  then  lowered  in  a  near 
vertical  mode  through  the  water  column.  If  a  ship  with  lifting  capabil¬ 
ities  were  available,  near-vertical  orientation  of  a  current  meter  coula 
be  obtained  by  attaching  a  large  weight  to  it's  base. 
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As  an  alternative  to  the  Savonius  rotor,  current  measurements  might 
be  obtained  through  other  means.  In  tnis  regard,  one  device  which  has 
be3n  recently  explored  and  shows  promise  for  measuring  flows  in  estuaries 
is  an  acoustic  current  meter.  Because  such  a  device  has  no  exposed 
moving  parts,  the  proolem  of  short-term  fouling  is  eliminated. 

If  additional  tests  were  to  be  performed  with  the  present  spar,  it  is 
suggested  that  they  be  made  in  deeper  water.  Since  the  portion  of  the 
spar  which  is  above  the  water  surface  serves  only  as  a  daymark,  or  plat¬ 
form  from  which  to  mount  one,  the  present  spar  might  be  modified  so  that 
a  small  daymark  shape  could  be  attached  to  the  top  using  small  diameter 
aluminum  tubing  such  as  that  used  for  airplane  struts.  In  this  way 
nearly  all  of  the  spar's  length  would  be  submerged  and  provide  a  large 
buoyant  moment.  As  a  result  much  of  the  large  degree  of  sensitivity  that 
this  spar  showed  to  small  depth  changes  would  be  eliminated. 
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User  Manual  for  the  Computer  Model,  ASB1 


< 


USERS  MANUAL  FOR  ThE  COMPUTER  MODEL,  AS81 

Introduction 

ft  program  has  been  developed  which  computes  the  total  force,  resul¬ 
tant  moment  and  if  articulated,  the  list  angle  on  a  spar  device  being 

acted  upon  by  wind,  current  and  wave  forces.  Consisting  of  a  main  pro¬ 
gram.  MAIN,  and  19  internally  calleo  subroutines,  this  program  written  in 
Fortran  IV  has  been  run  successfully  on  the  University  of  Rhode  Island’s 
ITELrtS  5  digital  computer.  Every  effort  has  been  made  to  make  this  pro¬ 
gram  as  general  as  possible.  Consequently,  a  wide  variety  of  input  data 
as  outlined  below  may  be  used.  Because  of  the  large  number  of  arrays 
which  may  be  used  during  computation,  it  is  necessary  to  allocate  320k 

oytes  of  core  to  execute  this  program.  Unfortunately,  as  a  result  of  tne 

wide  range  of  possible  input  data  as  well  as  tne  iterative  nature  of  tne 

program,  it  is  not  possible  to  estimate  the  amount  of  time  required  to 
compile  and  execute  this  analytical  model. 

This  orogram  may  be  thought  of  as  having  two  distinct  phases;  a  con¬ 
trol  phase  and  a  computational  phase.  Program  MAIN  acts  in  a  control 

capacity,  reading  input  data,  checking  its  validity  and  calling  the 

appropriate  suo routines.  During  the  computational  procedure,  program 
MAIN  keeps  track  of  tne  internal  counters  present  in  the  orogram  and 
ensures  that  the  iterative  technique  used  to  determine  the  spar  list 
angle  converges.  Once  the  appropriate  force,  moment  and  list  angle 

information  has  been  computed  for  a  given  case,  program  main,  checxs  to 
see  if  there  are  more  cases  and  if  so  resets  all  counters  to  tneir 

initial  value. 

The  computational  phase  of  the  program  consists  of  the  subroutines. 
All  of  the  actual  calculations  take  place  in  these  subroutines,  tnere 
being  at  least  one  subroutine  for  each  of  tne  major  steps  shown  in  the 


flowchart  (Figure  II-3).  These  computational  operations  are  repeated  in 
a  cyclical  manner  until  the  convergence  criteria  specifiea  in  the  main 
program,  MAIN,  is  satisfied. 


Input  Data 

The  English  system  of  units  (pounds-mass,  feet,  seconds)  is  assumed 
for  all  input  data  unless  otherwise  specified.  Since  most  of  the 
environmental  ana  structural  data  available  to  users  of  this  program  is 
described  using  this  system  of  units,  it  was  chosen  to  maKe  running  this 
program  as  easy  as  possible.  All  of  the  input  data  is  read  by  the  main 
program,  MAIN,  where  it  is  checked  to  ensure  that  each  value  is  within 
its  allowable  range.  Should  an  input  value  be  specified  outside  of  this 
range,  a  message  to  this  effect  is  Printed  out  and  the  Program  is  termi¬ 
nated.  Input  data  may  be  grouped  into  three  catagories;  general  informa¬ 
tion,  soar  structural  information,  and  force  information  3nd  should  be 
read  into  the  program  in  the  following  order.  Unless  otherwise  speci¬ 
fied,  an  F10.0  format  is  assumed  for  all  input  variaoles. 

Program:  MAIN 

Input  Variables 

General  Information 

NCASE  Numoer  of  cases  contained  in  the  inout  data  set  (Integer  *3) 

I  WAT  Water  type.  This  value  is  used  to  define  the  fluid  density 

and  viscosity.  It  should  be  set  equal  to  1  if  the  buoy  is 
to  be  located  in  salt  water  or  0  if  this  is  not  the  case. 
(Integer  *  1) 


2 


OEPTH 


Still  water  depth. 


Spar  Structural  Information 

SPL.SPW.SPD  Spar  length;  soar  width;  spar  depth. 

5M,SMO,SWT  Spar  material  (Alpna-numeric  *  10  J ;  spar  material  density; 

spar  wall  thicKness.  If  the  spar  is  a  solid  structure, 
then  tne  spar  wall  thickness  value  should  oe  set  equal  to 
0.00. 

IVANE  Orienting  vane  flag.  This  value  should  be  set  equal  to  1 

if  an  orienting  vane  is  present,  otherwise  this  value 
should  be  set  equal  to  0  and  the  next  three  vane  informa¬ 
tion  data  cards  omitted.  (Integer  *  1). 

VL ,  VHT , VT  Vane  length;  vane  height;  vane  thickness. 

VM ,  VMD  Vane  material  (Alpha-numeric  *  10 ) ,  vane  material  density. 

VHGT  Vane  height.  This  is  the  height  from  the  bottom  of  tne 

spar  to  the  bottom  sage  of  the  orienting  vane. 

IART  Articulation  flag.  If  the  buoy  is  articulated  aocut  its 

base  then  this  value  should  be  set  equal  to  0.  If  this 
value  is  not  equal  to  1,  then  it  is  assumed  that  the  spar 
acts  as  a  pile  structure  which  is  aligned  perpendicular  to 
the  bottom  and,  therefore,  no  list  angles  or  rignting 
moments  are  computed.  (Integer  *  1). 

ANCHXjANCHY  Coordinates  of  the  anchor  attachment  point  relative  to  tne 
center  of  the  bottom  face  of  the  spar. 

Force  Information 

NWOP  lumber  of  wind  data  points.  This  value  establishes  the 

number  of  elemental  areas  into  wnich  that  portion  of  the 
spar  which  is  above  the  still  water  surface  is  divided. 


WSPO,WGIR 


WLEEV, FETCH 

NCOP 


I  CUR 


UCVEL.UCDIR 


Because  of  dimension  limitations  on  the  wind  arrays,  this 
value  should  always  be  less  than  100.  if  tnis  value  is  set 
equal  to  0,  then  the  next  wind  information  data  card  must 
oe  omitted.  (Integer  *3) 

Wind  speed  (in  miles  per  hour);  wind  direction  (in  degrees 
magnetic) 

elevation  at  which  the  wind  speed  was  measured;  fetch 
length. 

Number  of  current  data  points.  This  value  establishes  the 
number  of  elemental  areas  into  wnich  the  suomerged  portion 
of  the  spar  is  divided  for  the  purpose  of  computing  current 
velocity  information.  This  value  should  not  exceed  500 
because  of  dimension  limitations  on  tne  current  arrays.  If 
this  value  is  set  equal  to  0,  it  is  assumed  that  no  current 
forces  are  present  ana  thus  the  next  tnree  current  informa¬ 
tion  data  cards  must  be  omitted.  (Integer  *  3) 

Uniform  current  flag.  This  value  determines  wnether  a 
single  set  of  variables  or  an  array  of  current  information 
are  the  next  input  avariables.  This  value  indicates 
whether  the  current  velocity  profile  is  uniform  or  not.  if 
this  value  is  set  to  1  then  the  following  data  card  snould 
contain  the  variaDles  UCVEL  and  UCOiR.  If  this  value  is  0 
then  the  subsequent  data  cards  should  contain  the  current 
velocity  Profile  VEL,  DIR,  OEP.  (Integer  *  1) 

Uniform  current  velocity;  uniform  current  direction.  Tnese 


arrav  elements  are  read  in  only  when  tne  L-UR  value  has 
been  set  equal  to  1. 


VEl,0IR,0£P 


NWAOP 


IWAVE 


Current  velocity;  current  direction;  current  depth.  These 
array  values  are  read  in  when  the  1CUR  variable  has  oeen 
set  equal  to  a  value  other  than  1.  having  a  length  equal 
to  NCOP  +  1,  each  row  of  this  array  contains  the  current 
speed  and  direction  at  an  associated  depth.  Beginning  with 
a  water  depth  of  O.QQ,  this  current  profile  information 
must  oe  supplied  at  equally  spaced  depth  increments. 

Number  of  wave  data  points.  This  value  establishes  the 
number  of  elemental  areas  into  which  the  submerged  portion 
of  the  spar  (as  measured  from  the  still  water  surface)  is 
divided  for  the  purpose  of  computing  wave  particle  velocity 
and  acceleration  information.  Because  of  dimension  limita¬ 
tions  on  the  wave  arrays,  this  value  cannot  exceed  500. 

Furth,r,  if  both  wave  and  currene  forces  are  present  in  an 
analysis,  this  value  must  be  less  than  or  equal  to  NCOP. 
If  the  nwagp  value  specified  is  non-zero  and  less  than  NCuP 
then  it  is  increased  to  equal  NCOP  and  a  message  to  the 
effect  is  printed,  if  NWAOP  is  set  to  0,  no  wave  forces 

are  present  and  the  next  four  wave  data  cards  must  be 

omitted.  (Integer  *  1) 

Wave  flag.  If  the  significant  and  average  wave  heights  ana 
wave  periods  are  to  be  supplied,  this  value  should  be  set 
to  1.  If  this  value  is  other  than  1  then  this  wave  infor¬ 
mation  is  internally  calculated  in  subroutine  WAVSMB  and, 
therefore,  the  next  wave  data  card  should  oe  omitted. 

(Integer  *  1) 


TSIG,TAVG,  Significant  wave  period;  average  wave  period; 


HGTSIG.HGTAVG  significant  wave  height;  average  wave  height.  These  values 


I  TIME 


ITHERY 


ION 


Additional 

NEWPO 


are  read  in  only  when  the  IWAVE  parameter  is  equal  to  1. 

Time  step  increment.  This  value  establisnes  the  numoer  of 
time  increments  into  wnich  the  average  wave  period  is 
divided.  The  total  force  on  the  spar  and  its  resulting 
list  angle  is  computed  at  each  of  these  time  steps. 
(Integer  *  3) 

Wave  theory  flag.  This  value  should  be  set  equal  to  1  if 
linear,  (Airy)  wave  theory  is  to  be  used  to  compute  the 
wave  particle  motions  and  surface  profile.  A  value  of  2  is 
used  to  specify  that  Stokes  second  order  wave  theory  is  to 
be  used  to  compute  these  wave  parameters.  (Integer  *  1) 

Integer  output  number.  This  value  should  be  set  equal  to  1 
if  a  listing  of  all  of  the  intermediate  calculations  used 
to  compute  the  resultant  list  angle  is  desired.  if  this 

value  is  other  than  1,  no  such  listing  will  be  produced. 
(Integer  *  1) 

General  Information 

New  physical  data.  This  value  is  only  read  in  wnen  the 
number  of  cases  in  the  input  data  set  is  greater  than  1. 
If  new  general  and  soar  structural  information  is  to  be 
forthcoming  tnen  this  value  should  be  set  to  1,  otherwise 
it  is  assumed  that  only  the  force  information  is  to  be 
varied  and  the  value  read  in  is  assigned  to  the  variable 


NWGP. 


THE  SUbRGUT iNES 


The  suoroutines  described  below  are  listed  in  the  order  in  which  they 
are  called  by  the  program,  MAIN.  The  first  four  subroutines  are  C3llea 
only  once  by  the  main  program  because  the  calculations  involvec  are 
independent  of  the  spar's  list  angles. 

Subroutine  WEIGHT 

The  WEIGHT  subroutine  computes  the  weight  of  the  entire  spar  cevice 
in  both  air  3nd  water.  If  tne  spar  has  been  previously  identified  oy  tne 
1  T  parameter  as  being  articulated  then  its  total  buoyancy  is  deter¬ 
mined.  3ased  upon  this  information,  a  check  is  tnen  made  to  ensure  that 
the  specified  spar  design  acting  as  a  free  body  in  the  absence  of  any 
environmental  forces  will  float.  Should  the  dimensions  of  the  spar  have 
been  specified  so  that  a  particular  design  will  noc  float,  a  message  to 
this  effect  is  printed  and  the  program  terminated. 

Subroutine  WNOSPO 

The  WNOSPO  subroutine  computes  the  wind  soeed  profile  from  the  wind 
velocity  input  data  based  uoon  the  Pranctl-von  Karman  universal  distribu¬ 
tion  law. 

Subroutine  WAV SMB 

The  significant  wave  height  and  significant  wave  period  of  a  wave  is 
computed  in  subroutine  WAVSM6  using  the  S.M.6.  method.  Assuming  that  the 
wave  height  and  the  square  of  the  wave  period  are  both  Rayleigh  distrio- 
uted,  average  values  for  these  two  parameters  are  also  determined. 
(Ippen,  1965).  These  average  values  are  subsequently  used  in  the  wave 
force  calculations. 


Suoroutine  WAVE 


Subroutine  WAVE  is  called  by  the  program,  MAIN,  to  compute  oasic 
information  about  the  waves  acting  on  the  spar  including  the  wavelength, 
wave  celerity,  wave  frequency  and  the  wave  number.  A  wave  amplitude 
value  based  upon  the  previously  selected  wave  theory  is  also  computed. 

Subroutine  AREAS 

In  this  analysis,  the  total  area  exposed  to  wind  forces  is  determined 
by  the  distance  from  mean  sea  level  to  the  top  of  the  spar;  the  total 
current  area  is  determined  by  the  distance  from  mean  sea  level  to  the 
oottom  of  the  spar  and  the  total  area  exposed  to  wave  particle  forces  is 
determined  by  tne  distance  from  the  free  surface  to  the  bottom  of  tne 
spar.  These  distances,  which  are  computed  in  subroutine  AREAS,  vary  as  a 
function  of  the  spar's  list  angle.  Further,  since  the  distribution  of 
the  driving  force  velocities  and  accelerations  may  nave  large  spatial 
variations  in  the  vertical,  subroutine  AREAS  divided  tnese  tnree  areas 
into  many  equal  elemental  areas.  The  exact  number  of  elemental  areas 
into  which  an  area  exposed  to  a  driving  force  is  divided  is  equal  to  one 
less  than  the  number  of  data  points  assigned  to  that  force  (i.e.,  one 
less  than  NWGP,  nco°Anq,  nwaDP). 


Subroutine  C0EAR 

8ecause  the  environmental  forces  in  this  model  are  defined  with 
respect  to  a  fixed  reference  point  (i.e.,  magnetic  north),  it  is  neces¬ 
sary  to  define  tne  orientation  of  the  spar  in  relationship  to  this  point 
in  order  that  the  relative  forces  acting  on  the  spar  be  established. 
Subroutine  CBEAR  is  called  to  compute  this  orientation  when  the  buoy 
configuration  under  study  has  a  splitter  plate  or  other  vane.  This 


subroutine  assumes  that  a  buoy's  orientation  is  defined  by  the  current 
direction  at  the  centroid  of  its  vane. 


The  Velocity  Subroutines 

Once  the  orientation  of  the  spar  has  been  established  with  respect  to 
magnetic  north,  the  driving  force  velocities  normal  to  the  spar  in  the  X 
and  Y  directions  are  determined. 

Subroutine  RWVEL 

Subroutine  RWVEL  calculates  the  wind  velocities  normal  to  the  spar  in 
the  x  and  Y  directions.  In  this  subroutine,  velocities  are  computed  at 
the  center  of  each  elemental  wind  area. 


Subroutine  RCVEL 

Current  \elocities  normal  to  the  spar  in  the  X  and  Y  directions  are 
ccmouted  in  subroutine  RCVEL.  These  velocities  are  comDuted  at  the  too 
and  bottom  of  each  elemental  current  area  and  then  averaged  linearly  to 
determine  the  velocity  at  the  center  of  each  area. 

Subroutine  RWAVA 

In  suorouine  RWAVA  the  water  particle  velocities  and  accelerations 
normal  to  the  spar  in  the  x  and  Y  directions  are  computed.  These  values 
are  computed  at  the  top  and  bottom  of  each  elemental  wave  particle  area 
and  then  averaged  linearly  to  determine  the  velocity  and  acceleration 
terms  at  the  center  of  each  area. 
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Subroutine  WPV 


In  oraer  to  properly  compute  the  total  force  due  to  current  and  waves 
acting  on  the  spar,  it  is  necessary  to  comoute  the  total  water  particle 
velocity  by  summing  these  individual  components  at  each  depth  of  interest. 

"When  current  is  present,  fluid 
forces  may  still  be  evaluated  accord¬ 
ing  to  Morrison's  formula  (Myers, 

1969).  This  is  acnieved  by  consider¬ 
ing  fluid  particle  velocity  as  the 
vector  sum  current  velocity  and  wave 
induced  particle  velocity.  Since  the 
drag  force  is  non-linear,  it  cannot  be 
regarded  as  a  simple  superposition  of 
current  and  wave  drag  forces."  (Wu 
and  Tung,  1975) 

Subroutine  WPV  combines  these  current  and  wave  particle  velocities  to 
determine  the  water  particle  velocity  and  associated  water  particle  area 
at  each  depth  of  interest. 


The  Orag  Coefficient  Subroutines 

Subroutine  DRAGC 

Subroutine  ORAGC  is  called  by  program  MAIN  when  the  spar  being 
studied  has  a  cylindrical  cross-section  vi.e.,  SPW  =  :>PD).  In  this 
subroutine,  drag  coefficients  are  computed  at  each  elemental  area  from 
Reynolds  number  determinations  according  to  the  U.5.  Army  Cores  of 
Engineers  Recommended  design  Curve.  (U.S.  Army  Corps  of  Engineers,  1973). 

Subroutine  ORAGX 

If  the  soar  under  study  does  not  have  a  cylindrical  cross-section 
then  drag  coefficients  for  driving  force  velocities  in  the  x-direction  to 
each  elemental  area  are  calculated  in  subroutine  ORAGX  pased  upon  empiri¬ 
cally  derived  relationships.  (Hoemer,  1965) 


Subroutine  DRAGY 


No  empirical  relationship  describing  the  drag  coefficient  on  an 
elliptical  body  whose  major  axis  was  normal  to  the  flow  could  be  found  in 
the  literature.  Consequently,  this  subroutine  assumes  that  such  a  snape 
may  be  considered  a  flat  plate  and,  therefore,  assigns  a  drag  coefficient 
value  of  2.0  to  driving  force  components  acting  in  this  airection. 

Subroutine  FORCE 

The  wind  and  water  particle  normal  forces  on  each  of  these  elemental 
areas  is  computed  in  subroutine  FORCE  using  the  Morison  equation  and 
velocity,  acceleration,  area  and  drag  coefficient  information  from 
previous  subroutines.  Subsequently,  the  total  wind  and  water  particle 
normal  force  on  the  spar  is  determined  by  summing  each  set  of  these 
elemental  forces. 

Subroutine  QTMGM 

Sunroutine  OTMOM  computes  the  total  wind  and  water  particle  overturn¬ 
ing  moments  acting  on  the  spar  in  the  X  and  Y  directions.  These  values 
are  subsequently  added  togetner  so  that  the  total  overturning  in  each  of 
these  directions  is  known. 

Subroutine  RMQM 

This  subroutine  calculates  the  total  righting  moment  acting  on  a  spar 


in  the  x  and  Y  directions  based  upon  buoyancy  considerations  and  the 
total  soar  weight. 


The  Intermediate  Outout  Subroutines 

For  tne  purpose  of  checking  the  final  output  or  the  computations 
performed  in  a  subroutine,  the  intermediate  output  subroutines  when 
called  by  program  MAIN  print  the  results  of  all  significant  calculations 
used  to  determine  the  total  force,  resulting  moment  and  list  of  angle  of 
a  spar. 

Subroutine  INOUT  1 

This  subroutine  prints  the  results  of  all  the  significant  calcula¬ 
tions  performed  in  the  following  subroutines:  AREAS,  C3EAR,  RWVEl, 
RCVEL,  RWAVA,  DRAGC ,  DRAGX,  and  DRAGY. 

Suo routine  INOUT  2 

This  subroutine  prints  the  results  of  all  the  significant  calcula¬ 
tions  performed  in  the  following  subroutines:  FORCE,  OTMGM,  AND  RMOM. 

Subroutine  anchor 

Subroutine  ANCHOR  determines  the  maximum  buoyant  force  acting  on  the 
soar  during  one  wave  period.  Subsequently,  the  minimum  standardized 
anchor  size  needed  to  resist  this  buoyant  force  in  determined  and  tne 
factor  of  safety  associated  with  this  anchor  selection  is  computed. 


APPENDIX  C 


Oesign  Orawings  of  the  Field  Test 


